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Precision electroweak data presently favors a weakly-coupled Higgs sector as the mechanism 
responsible for electroweak symmetry breaking. Low-energy supersymmetry provides a natural 
Cn| ■ framework for weakly-coupled elementary scalars. In this review, we summarize the theoretical 

properties of the Standard Model (SM) Higgs boson and the Higgs sector of the minimal super- 
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\ symmetric extension of the Standard Model (MSSM). We then survey the phenomenology of the 

SM and MSSM Higgs bosons at the Tevatron, LHC and a future e + e~ linear collider. We focus on 
Qh. the Higgs discovery potential of present and future colliders and stress the importance of precision 

r- 1 \ measurements of Higgs boson properties. 

> : 

X ; 1 Introduction — Origin of Electroweak Symmetry Breaking 

Deciphering the mechanism that breaks the electroweak symmetry and generates the masses of the 
known fundamental particles is one of the central challenges of particle physics. The Higgs mechanism [jlj 
in its most general form can be used to explain the observed masses of the and Z bosons as a 
consequence of three Goldstone bosons (G ± and G°) that end up as the longitudinal components of 
the gauge bosons. These Goldstone bosons are generated by the underlying dynamics responsible for 
electroweak symmetry breaking. However, the fundamental nature of this dynamics is still unknown. 
Two broad classes of electroweak symmetry breaking mechanisms have been pursued theoretically. In 
one class of theories, electroweak symmetry breaking dynamics is weakly-coupled, while in the second 
class of theories the dynamics is strongly-coupled. 

The electroweak symmetry breaking dynamics that is employed by the Standard Model posits a 
self- interacting complex doublet of scalar fields, which consists of four real degrees of freedom ||. 
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Renormalizable interactions are arranged in such a way that the neutral component of the scalar doublet 
acquires a vacuum expectation value, v = 246 GeV, which sets the scale of electroweak symmetry 
breaking. Consequently, three massless Goldstone bosons are generated, while the fourth scalar degree 
of freedom that remains in the physical spectrum is the CP-even neutral Higgs boson (hsu) of the 
Standard Model. It is further assumed in the Standard Model that the scalar doublet also couples 
to fermions through Yukawa interactions. After electroweak symmetry breaking, these interactions are 
responsible for the generation of quark and charged lepton masses. This approach is an example of weak 
electroweak symmetry breaking. Assuming that rrih SM < 200 GeV, all fields remain weakly interacting 
at energies up to the Planck scale. In the weakly-coupled approach to electroweak symmetry breaking, 
the Standard Model is very likely embedded in a supersymmetric theory |3] in order to stabilize the 
large gap between the electroweak and the Planck scales in a natural way @||. These theories predict 
a spectrum of Higgs scalars with the properties of the lightest Higgs scalar often resembling that of 
the Standard Model (SM) Higgs boson. 

Alternatively, strong breaking of electroweak symmetry is accomplished by new strong interactions 
near the TeV scale 0. More recently, so-called "little Higgs models" have been proposed in which 
the scale of the new strong interactions is pushed up above 10 TeV [§], and the lightest Higgs scalar 
resembles the weakly-coupled SM Higgs boson. In a more speculative direction, a new approach to 
electroweak symmetry breaking has been explored in which extra space dimensions beyond the usual 
3 + 1 dimensional spacetime are introduced |9| with characteristic sizes of order (TeV) -1 . In such 
scenarios, the mechanisms for electroweak symmetry breaking are inherently extra-dimensional, and 
the resulting phenomenology may be significantly different from the usual approaches mentioned above. 
All these alternative approaches lie outside the scope of this review. 

Although there is as yet no direct evidence for the nature of electroweak symmetry breaking dynam- 
ics, present data can be used to discriminate among the different approaches. For example, precision 
electroweak data, accumulated in the past decade at LEP, SLC, the Tevatron and elsewhere, strongly 
support the Standard Model with a weakly-coupled Higgs boson [Kj. Moreover, the contribution of 
new physics, which can enter through W ± and Z boson vacuum polarization corrections, is severely 
constrained. This fact has already served to rule out several models of strongly-coupled electroweak 
symmetry breaking dynamics. The Higgs boson contributes to the W ± and Z boson vacuum polariza- 
tion through loop effects, and so a global Standard Model fit to the electroweak data yields information 
about the Higgs mass. The results of the LEP Electroweak Working Group analysis shown in fig. |l](a) 
yield PTJfl : m hsM = 81 GeV, and provides a 95% CL upper limit of m hsM < 193 GeV. These re- 



sults reflect the logarithmic sensitivity to the Higgs mass via the virtual Higgs loop contributions to 
the various electroweak observables. The 95% CL upper limit is consistent with the direct searches at 



LEP [11] that show no conclusive evidence for the Higgs boson, and imply that rrih SM > 114.4 GeV at 



95% CL. Fig. |l](b) exhibits the most probable range of values for the SM Higgs mass [0. This mass 
range is consistent with a weakly-coupled Higgs scalar that is expected to emerge from the Standard 
Model scalar dynamics. 
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Figure 1: (a) The "blueband plot" shows A\ 2 
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as a function of the Standard Model Higgs mass [hoi . The 



solid line is a result of a global fit using all data; the band represents the theoretical error due to missing higher order 
corrections. The rectangular shaded region shows the 95% CL exclusion limit on the Higgs mass from direct searches 
at LEP jy]. (b) Probability distribution function for the Higgs boson mass, including all available direct and indirect 
data The probability is shown for 1 GeV bins. The shaded and unshaded regions each correspond to an integrated 
probability of 50% 



There are some loopholes that can be exploited to circumvent this conclusion. It is possible to 
construct models of new physics where the goodness of the global Standard Model fit to precision 
electroweak data is not compromised while the strong upper limit on the Higgs mass is relaxed. In 



particular, one can construct effective operators |13|,|T4| or specific models |15j of new physics where 
the Higgs mass is significantly larger, but the new physics contributions to the and Z vacuum 
polarizations, parameterized by the Peskin-Takeuchi \W\ parameters S and T, are still consistent with 
the experimental data. In addition, some have argued that the global Standard Model fit exhibits 
possible internal inconsistencies |T7|], which would suggest that systematic uncertainties have been 
underestimated and/or new physics beyond the Standard Model is required. Thus, although weakly- 
coupled electroweak symmetry breaking seems to be favored by the strong upper limit on the Higgs 
mass, one cannot definitively rule out all other approaches. 

Nevertheless, one additional piece of data is very suggestive. Within the supersymmetric extension 
of the Standard Model, grand unification of the electromagnetic, the weak and the strong gauge interac- 
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tions can be achieved in a consistent way, strongly supported by the prediction of the electroweak mixing 
angle at low energy scales with an accuracy at the percent level . The significance of this predic- 

tion is not easily matched by other approaches. For example, in strongly-coupled electroweak symmetry 
breaking models, unification of couplings is not addressed per se, whereas in extra-dimensional models 
it is often achieved by introducing new structures at intermediate energy scales. Unless one is willing to 
regard the apparent gauge coupling unification as a coincidence, it is tempting to conclude that weak 
electroweak symmetry breaking with low-energy supersymmetry is the preferred mechanism, leading 
to an expected mass of the lightest Higgs boson below 200 GeV (less than 135 GeV in the simplest 
super symmetric models), and a possible spectrum of additional neutral and charged Higgs bosons with 
masses up to of order 1 TeV. 

Henceforth, we shall assume that the dynamics of electroweak symmetry breaking is a result of 
a weakly-coupled scalar sector. The Standard Model is an effective field theory and provides a very 
good description of the physics of elementary particles and their interactions at an energy scale of 
(9(100) GeV and below. However, there must exist some energy scale, A, at which the Standard Model 
breaks down. That is, the Standard Model is no longer adequate for describing the theory above A, 
and degrees of freedom associated with new physics become relevant. In particular, we know that 
A < Mpl, since at an energy scale above the Planck scale, Mpl — 10 19 GeV, quantum gravitational 
effects become significant and the Standard Model must be replaced by a more fundamental theory that 
incorporates gravity. (Similar conclusions also apply to recently proposed extra-dimensional theories in 
which quantum gravitational effects can become significant at energies scales as low as (9(1 TeV) 0.) 
Of course, it is possible that new physics beyond the Standard Model exists at an energy scale between 
the electroweak and Planck scale, in which case the value of A might lie significantly below Mpl{] 

The value of the Higgs mass itself can provide an important constraint on the value of A. If rrih SM is 



too large, then the Higgs self-coupling blows up at some scale A below the Planck scale p(J . If ?rih SM is 
too small, then the Higgs potential develops a second (global) minimum at a large value of the scalar 
field of order A ||21|| . Thus new physics must enter at a scale A or below in order that the global minimum 



of the theory correspond to the observed SU(2)xU(l) broken vacuum with v = 246 GeV. Thus, given 
a value of A, one can compute the minimum and maximum Higgs mass allowed. The results of this 
computation (with shaded bands indicating the theoretical uncertainty of the result) are illustrated in 
fig. ^(a) P2| . Consequently, a Higgs mass range 130 GeV < m hsM < 180 GeV is consistent with an 
effective Standard Model that survives all the way to the Planck scale. f\ 

However, the survival of the Standard Model as an effective theory all the way up to the Planck 



^or example, the recent experimental evidence for neutrino masses of order 10 -2 eV or below cannot be strictly 
explained in the Standard Model. Yet, one can easily write down a dimcnsion-5 operator responsible for neutrino masses 
that is suppressed by v/A. If m„ ~ 10 -2 eV, then one obtains as a rough estimate A < 10 15 GeV. 

2 The constraint on A due to vacuum stability in fig. |^ is less stringent if one allows for the electroweak vacuum to 
be metastable, with a lifetime greater than the age of the universe. An analysis of rcf. Q finds that for a sufficiently 
long-lived electroweak vacuum, the Higgs mass lower limit of 130 GeV just quoted is reduced to about 115 GeV. 
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Figure 2: (a) The upper |^0| and the lower |H| Higgs mass bounds as a function of the energy scale A at which the 
Standard Model breaks down, assuming Mt = 175 GeV and a s (mz) = 0.118, taken from ref. p^j. The shaded areas above 
reflect the theoretical uncertainties in the calculations of the Higgs mass bounds, (b) Following ref. Jl4j], a reconsideration 
of the A vs. Higgs mass plot with a focus on A < 100 TeV. Precision electroweak measurements restrict the parameter 
space to lie below the dashed line, based on a 95% CL fit that allows for nonzero values of S and T and the existence of 
higher dimensional operators suppressed by v 2 /A 2 . The unshaded area has less than one part in ten fine-tuning. 



scale is unlikely based on the following "naturalness" [|J argument. In an effective field theory, particle 
masses and dimensionless couplings of the low-energy theory are calculable in terms of parameters of 
a more fundamental theory that describes physics at the energy scale A. All low-energy couplings and 
fermion masses are logarithmically sensitive to A. In contrast, scalar squared-masses are quadratically 
sensitive to A. Thus, in this framework, the observed Higgs mass (at one-loop) has the following form: 

2 / 2 \ ^9 ^ / -i \ 

™h SM = ( m h)0 + , (1) 

where (mh)o is a parameter of the fundamental theory, g is an electroweak coupling and k is a con- 
stant, presumably of 0(1), that is calculable within the low-energy effective theory. Because these 
two contributions arise from independent sources, it is very unlikely that the magnitude of is 
significantly smaller than either of the two terms. That is, the "natural" value for the physical scalar 
squared- mass is at least of order g 2 A 2 /16tt 2 . In order for this value to be consistent with the requirement 
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that the Higgs mass is of order the electroweak symmetry breaking scale (as required from unitarity 



constraints [24,2a]), the value of A must satisfy 



A ~ 4?rmfesM rsjQCl TeV). (2) 
9 

If A is significantly larger than 1 TeV (often called the hierarchy problem in the literature), then the 
only way to generate a Higgs mass of 0(rriz) is to have an "unnatural" cancellation between the two 
terms of eq. (p]). This seems highly unlikely given that the two terms of eq. (|]) have completely different 
origins. The requirement of A ~ 0(1 TeV) as a condition for the absence of fine-tuning of the Higgs 
mass parameter is nicely illustrated in fig. 0(b), taken from ref. 



A viable theoretical framework that incorporates weakly-coupled Higgs bosons and satisfies the 
constraint of eq. (Q) is that of "low-energy" or "weak-scale" supersymmetry J3|. If supersymmetry 
(which relates fermion and boson masses and interactions) is exact, then boson masses must exhibit the 
same logarithmically sensitivity to A as do the fermion masses. Since no supersymmetric partners of 
the Standard Model particles have been found, it follows that supersymmetry is not an exact symmetry 
of the fundamental particle interactions. Hence, in the framework of low-energy supersymmetry, A 
should be identified with the energy scale of supersymmetry-breaking. The naturalness constraint of 
eq. (§) is still relevant, which implies that the scale of supersymmetry breaking should not be much 
larger than a few TeV in order that the naturalness of scalar masses be preserved. Moreover, low-energy 
supersymmetry with a supersymmetry-breaking scale of (9(1 TeV) is precisely what is needed to explain 
the observed gauge coupling unification as previously noted. We conclude that a suitable replacement 
for the Standard Model is a supersymmetric extension of the Standard Model as the effective field theory 
of the TeV scale. One good feature of the supersymmetric approach is that the effective low-energy 
supersymmetric theory can be valid all the way up to the Planck scale, while still being natural! 

The physics of the Higgs bosons will be explored by experiments now underway at the upgraded 
proton-antiproton Tevatron collider at Fermilab and in the near future at the Large Hadron Collider 
(LHC) at CERN. Once evidence for electroweak symmetry breaking dynamics is obtained, a more com- 
plete understanding of the mechanism involved will require experimentation at a future e + e~ linear 
collider (LC) now under development. In this review we focus primarily on the theory and phenomenol- 
ogy of the Standard Model Higgs boson and the Higgs bosons of low-energy supersymmetry. In Section 0, 
we review the theoretical properties of the Standard Model Higgs boson, and exhibit its main branching 
ratio and production rates at hadron colliders and at the LC. The main Higgs boson search techniques 
at the Tevatron, LHC and the LC are described. In Section |3|, we examine the Higgs bosons of the min- 
imal supersymmetric Standard Model (MSSM). We summarize the tree-level properties of the MSSM 
Higgs sector and describe the most significant effects of the radiative corrections to the computation of 
the Higgs masses and couplings. We then exhibit the main branching ratios and production rates of the 
MSSM Higgs bosons and survey the phenomenology of the MSSM Higgs sector at the Tevatron, LHC 
and LC. A brief summary concludes this review in Section E|. 
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2 The Standard Model Higgs Boson 



In the Standard Model, the Higgs mass is given by: = \\v 2 , where A is the Higgs self-coupling 

parameter. Since A is unknown at present, the value of the Standard Model Higgs mass is not predicted. 
However, other theoretical considerations, discussed in Section [I], place constraints on the Higgs mass 
as exhibited in fig. ^|. In contrast, the Higgs couplings to fermions [bosons] are predicted by the theory 
to be proportional to the corresponding particle masses [squared- masses] . In particular, the SM Higgs 
boson is a CP-even scalar, and its couplings to gauge bosons, Higgs bosons and fermions are given by:0 

_m f _ 2m\ 2m 2 v 
9hff — — ) 9hvv — > 9hhvv — — 5— , \o) 

3 ^ m h SM 3 ^ m h S M 
ghhh = 2^ V = ~ ' Qhhhh — = > (4) 

where h = hsM, V = W or Z and v = 2m w /g = 246 GeV. In Higgs production and decay processes, the 
dominant mechanisms involve the coupling of the Higgs boson to the W ± , Z and/or the third generation 
quarks and leptons. Note that a h^gg coupling (g=gluon) is induced by virtue of a one-loop graph in 
which the Higgs boson couples to a virtual tt pair. Likewise, a /isM77 coupling is generated, although in 
this case the one-loop graph in which the Higgs boson couples to a virtual W + W~ pair is the dominant 
contribution. Further details of the SM Higgs boson properties are given in ref. 0. A review of the 
SM Higgs properties and its phenomenology, with an emphasis on the impact of loop corrections to the 
Higgs decay rates and cross-sections can be found in ref. p6| . 



2.1 Standard Model Higgs Boson Decay Modes 

The branching ratios for the main decay modes of a SM Higgs boson are shown as a function of Higgs 
boson mass in fig. ||] and f|(a), based on the results obtained using the HDECAY program [27]. For Higgs 



boson masses below 135 GeV, the decay h^m — * bb dominates, whereas above 135 GeV, the dominant 
decay mode is /ism WW^*' (below W + W~ threshold, one of the W bosons is virtual as indicated by 
the star). Above tt threshold, the branching ratio into top-quark pairs increases rapidly as a function of 
Higgs mass, reaching a maximum of about 20% at mh SM ~ 450 GeV. The total Higgs width is obtained 
by summing all the Higgs partial widths and is displayed as a function of Higgs mass in fig. Mb). 



The leading effects of the QCD corrections to the Higgs decay to quark pairs Eon can be taken 



into account by using the tree-level formula for the Higgs partial width (which depends on the quark 

3 The corresponding Feynman rules are obtained by multiplying the Higgs boson- VV couplings by ig^ and the other 
Higgs couplings of eqs. (||) and (Q) by a factor of —i. The appropriate combinatorial factors have been included. 

4 The formulae for the leading order QCD-corrections to T(/ism - * 9?) are nicely summarized in ref. [^9[| . The leading 
electroweak radiative corrections have been obtained in ref. [^p| . A useful summary of results, which includes new two- loop 
contributions to the radiatively corrected /ism — ► bb partial width, can be found in ref. |31j|. 
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Figure 3: Branching ratios of the dominant decay modes of the Standard Model Higgs boson as a function of Higgs mass 
for m/j SM < 200 GeV, taken from ref. |||. These results have been obtained with the program HDECAY [£7|, and include 
QCD corrections beyond the leading order |2^]. The shaded bands represent the variations due to the uncertainties in the 
input parameters: a s (Mf ) = 0.120±0.003, m b (M b ) = 4.22±0.05 GeV, m c (M c ) = 1.22±0.06 GeV, and M t = 174±5 GeV. 



mass), and identifying the quark mass with the running quark mass evaluated at the Higgs mass, 
mg(m/ lSM ). The running quark mass, mQ(mh SM ) is obtained from the MS mass, Wlq{Mq) [where Mq 
is the corresponding quark pole mass], by renormalization group evolution. The MS quark masses are 
obtained from fits to experimental data f5ci|] . Note that the large decrease in the charm quark mass 
due to QCD running is responsible for suppressing BR(cc) relative to BR(r + r~), in spite of the color 
enhancement of the former, thereby reversing the naively expected hierarchy. Below the corresponding 
two-body thresholds, the WW^*\ ZZ^ and t^t decay modes (where the asterisk indicates an off-shell 
particle) are still relevant as shown in fig. 

The h SM gg, /ism77 and /i SM Z7 vertices are generated at one-loop. The partial width for h SM — > gg is 
primarily of interest because it determines the gg — > h SM production cross-section. The /ism77 vertex is 
especially relevant both for the hswi — ► 77 discovery mode at the LHC and for the 77 — > hsM production 
mode at the LC operating as a 77 collider. 
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Figure 4: (a) Branching ratios of the Standard Model Higgs boson as a function of Higgs mass. Two-boson [fcrmion- 
antifermion] final states are exhibited by solid [dashed] lines. As compared with fig. a larger range of Higgs masses and 
branching ratios are shown, (b) The total width of the Standard Model Higgs boson is shown as a function of its mass. 
For comparison, we exhibit the widths of the two CP-even scalars, h and H of the MSSM for two different choices of 
MSSM parameters (tan (3 = 3 and 30 in the maximal mixing scenario; the onset of the H — > hh and H — ► tt thresholds 
in the tan/3 = 3 curve are clearly evident). The central values of a s , rnb(Mb) and m c (M c ) quoted in the caption of fig. || 
are employed in both (a) and (b). 



2.2 Standard Model Higgs Boson Production at Hadron Colliders 
2.2.1 Cross-sections at hadron colliders 



This section describes the most important Higgs production processes at the Tevatron (-^s 
and the LHC = 14 TeV). The relevant cross-sections are exhibited in figs. || and |6| f|32| 
Combining these Higgs production mechanisms with the decays discussed in Section BTT 
the most promising signatures. 

Due to the large luminosity of gluons at high energy hadron colliders, gg — > h SM is the Higgs 



2 TeV) 

Hi- 

one obtains 




production mechanism with the largest cross-section at the Tevatron and the LHC |37| , |38| . The two- 
loop, next-to-leading order (NLO) QCD corrections enhance the gluon fusion cross-section by about a 
factor of two [§8|,f^J] • The corresponding NLO differential cross-section (as a function of the Higgs boson 
Pt and rapidity) has also been obtained |4(|. Recently, the next-to-NLO (NNLO) QCD corrections have 
been evaluated ||41|| , and show a further enhancement of about 10% to 30% depending on the Higgs 
mass and center-of-mass energy of the collider. The remaining scale dependence and the effects of 
higher order terms not yet computed are estimated to give a theoretical uncertainty of 10-20%. The 
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dependence of the gluon fusion cross-section on different parton densities yields roughly an additional 
uncertainty of order 10%. 

The cross-section for qq — > W ± hsM (summed over both W charge states) is the second largest 
Higgs cross-section at the Tevatron for m hsM < 175 GeV. At the LHC, the W ± hsM cross-section is 
not as prominent over the Higgs mass range of interest. The corresponding qq — > ZhsM cross-section 
is roughly a factor of two lower than the corresponding W^hswi cross-section. The QCD corrections 
to o(qq — > Vhsyi) [V = W oi Z] coincide with those of the Drell-Yan process and increase the cross- 
sections by about 30% [42. 43|.|44|]. The theoretical uncertainty is estimated to be about 15% from the 



remaining scale dependence. The dependence on different sets of parton densities is rather weak and 
also leads to a variation of the production cross-sections by about 15%. 

Vector boson fusion is a shorthand notation for the full qq — > qqh^M process, where the quark and 
anti-quark both radiate virtual vector bosons (V*) which then annihilate to produce the Higgs boson. 
Vector boson fusion via ud — > duhsM (and its charge-conjugate process) is also possible. In figs. |] and [6|, 
all contributing processes are included, and the sum of all such contributions is labeled qq —> qqhsM 
for simplicity. The QCD corrections enhance the cross-section by about 10% E1,H5[. The vector boson 



fusion process is the second largest Higgs cross-section at the LHC; its cross-section approaches the 
99 hsM cross-section for rrih SM ~ 1 TeV. 

The cross-sections for gg, qq — > tihsM at the Tevatron and the LHC are displayed in figs. ^| and ||. 
The NLO QCD corrections to qq — > tihsM have recently been computed in refs. p6[ and ||47|| , and the 



corrections to gg — > ii/isM have been obtained in ref. |[47|| . Fig. |5] includes the complete NLO QCD 
corrections from ref. [|^| (at Tevatron energies, tt/isM production is dominated by the qq — > tihsM 
subprocess). The size of the QCD corrections depends sensitively on the choice of scale, /i, employed 
in the running coupling constant and parton distribution functions. Changes in /i can significantly 
modify the tree-level cross-sections, whereas the NLO-corrected cross-sections are rather insensitive to 
reasonable changes in /i. In fig. |5], the tree-level tihsM cross-section is shown for /i 2 = s (the square 
of the partonic center-of-mass energy), and is roughly a factor of two smaller than the corresponding 
cross-section with fi = m t . With respect to the latter choice, the NLO-corrections of ref. |^7| are rather 



small, typically of order 10-20% depending on the precise choice of 

The tree-level gg, qq — > bbh SM cross-section (as a function of m hsM ) shown in fig. ^| has been 
computed by fixing the scales of the parton distribution functions, the running coupling ot s and the 
running Higgs-bottom-quark Yukawa coupling (or equivalently, the running fe-quark mass) at the value 
of the corresponding /ism mass. In fig. |6|, the 6-quark pole mass is employed,^ while the scales for a s 
and the parton distribution functions were set equal to the partonic center-of-mass energy. 

By using the running b quark mass, one implicitly resums large logarithms associated with the QCD- 
corrected Yukawa coupling. Thus, the tree- level bbhsM cross-section displayed in fig. |5] implicitly includes 



5 The effect of using the 6-quark pole mass [Mf, ~ 5 GeV] as opposed to the running 6-quark mass ['7ib('rt/ lSM )] is to 
increase the cross-section by roughly a factor of two. 
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Figure 5: Higgs production cross-sections (in units of pb) at the Tevatron [y/s = 2 TeV], for the various production 
mechanisms as a function of the Higgs mass, taken primarily from refs. p3| and |3jJ. The full NLO QCD-corrected results 
are employed for the gluon fusion gg — > /ism, vector boson fusion qq — > V*V*qq — > hsMQQ (here, qq refers to both ud 
and qq scattering), Higgs-strahlung processes qq — > V* — » V/ism (where V = W^, Z), 66 — > /ism (taken from ref. p6[), 
and gg,qq — > hsMtt- Tree-level cross-sections are exhibited for gg,qq — ► hsyibb. In the latter case, the cross-section has 
been computed with a running 6-quark mass, a s and the parton distribution functions all evaluated at the corresponding 
Higgs mass. 



a part of the QCD corrections to the full inclusive cross-section. However, the most significant effect 
of the QCD corrections to bbh^M production arises from the kinematical region where the b quarks are 
emitted near the forward direction. In fact, large logarithms arising in this region spoil the convergence 
of the QCD perturbation series since a s ln(m^ gM /m^) ~ 0(1). These large logarithms (already present 
at lowest order) must be resummed to all orders, and this resummation is accomplished by the generation 
of the 6-quark distribution function [ji8l , fl9fl . Thus, the QCD-corrected fully inclusive bbhsM cross-section 



can be approximated by bb — >• hsM and its QCD corrections.^ The latter is also exhibited in fig. |5] and is 



6 This result, although correct in the far asymptotic regime (where *Js 3> TUh SM 3> mj), may still not be reliable for 
Higgs production at the Tevatron and LHC. In ref. |30), it is argued that even at the LHC for m^ SM = 500 GeV, sizable 
TOb-cffccts still remain and a(bb — > /ism) is an overestimate of the true QCD-corrected 66/ism cross-section. 
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Figure 6: Higgs production cross-sections (in units of pb) at the LHC [y/s = 14 TeV], for the various production 
mechanisms as a function of the Higgs mass, taken from rcf. ^4). The cross-section curves for gg, qq — ► /ism^ (which has 
not been updated to include the NLO calculation of ref. j47j) and gg, qq — ► hsM^b are based on a tree-level calculation with 
i-quark and 6-quark pole masses and a s and the parton distribution functions evaluated at the partonic center-of-mass 
energy. 



seen to be roughly an order of magnitude larger than the tree-level bbh^M cross-section [36]. Of course 



this result is not very relevant for the searches at hadron colliders in which transverse momentum cuts 
on the 6-jets are employed. Ultimately, one needs the QCD-corrected differential cross-section for bbh^M 
(as a function of the final state 6-quark transverse momentum) in order to do realistic simulations of the 
Higgs signal in this channel. However, if only one 6-quark jet is tagged, it may be sufficient to consider 



the process bg — > bhsM- The cross-section for bg — > fe/isM at lowest order can be found in ref. pi ]; the 
NLO QCD-corrected cross-section has been recently obtained in ref. ||52|| . For example, assuming that 
Pt > 15 GeV, and the pseudorapidity \rj\ < 2 for the observed fe-quark jet, the NLO cross-section at 
the Tevatron ranges from about 6 fb to 0.25 fb for 100 GeV < rnh SM < 200 GeV. Increasing the cuts to 
Pt > 30 GeV and \rj\ < 2.5 at the LHC yields a range of cross-sections from about 200 fb to 1.2 fb for 
100 GeV < m hsM < 500 GeV. 

Not shown in figs. |5] and || is the cross-section for inclusive double Higgs production (/ism^-sm + X). 
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Double Higgs production is not observable at the Tevatron, but may be possible to detect at the LHC 
given sufficient luminosity. The main contributions to double Higgs production in order of importance 
are: (i) gg -> h SM h SM ; (ii) VV -> h SM h SM ] and (iii) qq -> Vh SM hsM, where V = W or Z. The 
gluon-gluon fusion cross-section dominates by at least an order of magnitude, so we focus on this 
subprocess |53T J54|J^|56| . Including NLO QCD corrections, typical cross-sections for pp — > /t-sm^sm + X 
at y/s 



14 TeV range from about 40 fb to 10 fb for 100 < m hsM < 200 GeV |^J56 

h* 

1 4c 



There are two 
-> h SM h SM (via 

t*t* — > ZlSM^SM 



classes of diagrams that contribute to h SM h SM production via gluon fusion: gg n, SM 
the top-quark triangle diagram) which is sensitive to the triple-Higgs vertex, and gg — 
(via the top-quark box diagram) which is independent of the triple-Higgs vertex. Due to the relatively 
low cross-sections, it will be very challenging to extract information on the Higgs self-coupling parameter 
from LHC data. 

Finally, assuming that very forward protons can be tagged, diffractive production of Higgs bosons 
may provide a viable signal at the LHC (the corresponding cross-sections at the Tevatron are prob- 
ably too small) [0)0 -Q Such events are characterized by rapidity gaps, i.e., the absence of particle 
production between the forward protons and the centrally produced Higgs boson. In particular, for 
m hsM ~ 130 GeV, the exclusive process pp — > pph SM (with h SM — > bb) can be used at the LHC to obtain 
the Higgs mass with an accuracy of about 1% by measuring the invariant mass recoiling against the 



final state protons [59 



2.2.2 Standard Model Higgs Boson Searches at the Tevatron 

In the mass region of interest to the Tevatron Higgs search (100 GeV< m hsM < 200 GeV), the SM 
Higgs boson is produced most copiously via gg fusion, with a cross-section from about 1.0 — 0.1 pb. 
For m hsM < 135 GeV, the Higgs boson decays dominantly to bb. Since the cross-section for the QCD 
production of bb dijet events is orders of magnitude larger than the Higgs production cross-section, the 
gg — > hsM — ► bb channel is not a promising channel. For rrih SM > 135 GeV, the Higgs boson decays 
dominantly to H / H /< '*' ) (where W* is a virtual W), and the channel gg — > hsM — > WW^ is accessible 
to the Tevatron Higgs search [Q. 

Given sufficient luminosity, the most promising SM Higgs discovery mechanism at the Tevatron 
for rrih SM < 135 GeV consists of qq annihilation into a virtual V* (V = W or Z), where V* — > Vhsu 
followed by a leptonic decay of the V and /ism — ► bb |6l[ . The sum of the WhsM and ZhsM cross-sections 



is about 0.2 — 0.5 pb in the mass region of interest (100 < m hsM < 135 GeV), in which the dominant 
Higgs decay is /igM ~^ bb. These processes lead to three main final states, iubb, vvbb and £ + £~bb, that 
exhibit distinctive signatures on which the experiments can trigger (high px leptons and/or missing 
Et). The backgrounds are manageable and are typically dominated by vector-boson pair production, tt 
production and QCD dijet production. The signal efficiencies and backgrounds have all been estimated 



7 A critical comparison of various theoretical approaches to diffractive Higgs production can be found in ref. 
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with both the CDF Run 1 detector simulation and with the simple SHW simulation [02]. Monte Carlo 
estimates have been used for the backgrounds everywhere except in the vvbb channel, where there is 
a significant contribution from QCD bb dijet production. To be conservative, in ref. |62| the unknown 
QCD bb dijet background to the vvbb channel has been taken to be equal in size to the sum of all other 
contributing background processes. In addition, the separation of signal from background was optimized 
using neural network techniques [63], resulting in a demonstrable gain in the significance of the Higgs 
signal for the Ivbb and vvbb channels ||62||. The fe-tagging efficiencies and the bb mass resolution play 



a key role in determining the ultimate efficiency and background rejection. Much work remains, using 
real data studies, to optimize the performance in both these areas. 

For larger Higgs masses (rrih SM ^ 135 GeV) it is possible to exploit the distinct signatures present 
when the Higgs boson decay branching ratio to WW(*) becomes appreciable. In this case, there are final 
states with WW (from the gluon-fusion production of a single Higgs boson), and WWW and ZWW 
arising from associated vector boson-Higgs boson production. Three search channels were identified in 
ref. [B2| as potentially sensitive at these high Higgs masses: like-sign dilepton plus jets (•^ ± ^ ± jj) events, 
high-pr lepton pairs plus missing Et (£ + £~vv), and trilepton events. Of these, the first two 

were found to be most sensitive. The strong angular correlations of the final state leptons resulting 
from WW* is one of the crucial ingredients for these discovery channels |p^j65| , |60 l. 

The integrated luminosity required per Tevatron experiment as a function of Higgs mass to either 
exclude the SM Higgs boson at 95% CL or discover it at the 3cr or 5a level of significance, for the SHW 
analyses with neural net selection (see ref. |62| for details), is shown in Fig. 0(a). These results are based 
on the combined statistical power of both the CDF and D0 experiments. The bands extend from the 
neural net result on the low side upward in required integrated luminosity by 30% to the high side, as 
an indication of the range of uncertainty in the 6-tagging efficiency, bb mass resolution and background 
uncertainties. As the plots show, the required integrated luminosity increases rapidly with Higgs mass 



to 140 GeV, beyond which the high-mass channels play the dominant role. If rrih s 
lies just above the 95% CL exclusion limit achieved by LEP JTI]|, then 5 fb _1 



115 GeV, which 



of integrated luminosity 

per experiment would provide sufficient data to see a 3a excess above background. With 15 fb^ 1 of 
integrated luminosity per experiment, a 5a discovery of the Higgs boson would be possible. 

The final result shows that for an integrated luminosity of 10 fb _1 , if the SM Higgs boson mass lies 
beyond the discovery reach of the Tevatron, then one can attain a 95% CL exclusion for masses up to 
about 180 GeV. Moreover, if the SM Higgs happens to be sufficiently light (rrih SM < 125 GeV), then a 
tantalizing 3a effect will be visible with the same integrated luminosity. With about 25 fb _1 of data, 3a 
evidence for the Higgs boson can be obtained for the entire Higgs mass range up to 180 GeV. However, 
the discovery reach is considerably more limited for a 5a Higgs boson signal. With 30 fb _1 integrated 
luminosity delivered per detector, a 5a Higgs boson discovery may be possible for Higgs masses up to 
about 130 GeV, a significant extension of the LEP Standard Model Higgs search. The latter figure of 
merit is particularly significant when applied to the search for the lightest Higgs bosons of the MSSM. 



We address this case in Section 3.5.3 
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Figure 7: (a) The integrated luminosity required per Tevatron experiment, to either exclude a SM Higgs boson at 95% 
CL or observe it at the 3er or 5cr level, as a function of the Higgs mass [ |62| . (b) Expected 5<r discovery luminosity 
requirements for the SM Higgs boson at the LHC for one experiment, based on a study performed with CMS fast detector 
simulation, assuming statistical errors only ]7^] . The gg and W + W~ fusion processes are indicated respectively by the 
solid and dotted lines. 



Other Higgs signatures could help improve the sensitivity of the Higgs search at the Tevatron. In 
ref. ||62|| , channels containing the h^M — > t + t~ decay mode have not been studied, as the small branching 
ratio (less than 8%) makes the corresponding signal rates small. Still, a significant improvement of r- 
lepton identification could lead to a viable Higgs signal in the Higgs mass region 120 GeV< m hsM < 
140 GeV ||66| . Another possibility which has been explored is the detection of the Higgs boson via tihsM 



production (the Higgs boson is radiated off the top-quark), followed by h$M — ► bb. Initial studies j6T | 
suggested that this channel could be observable at the upgraded Tevatron for m hsM < 140 GeV, with a 
statistical significance comparable to the Higgs signals in the Wh^M and Zh^M channels. 

If a Higgs boson is discovered at the Tevatron, one can begin to measure some of its properties. The 
Higgs mass can be measured with an accuracy of about 2 GeV |68| |. However, the determination of the 



Higgs couplings to W and Z bosons and to bb will be mo del- dependent and rather crude. To improve 
and expand the possible Higgs measurements and determine its phenomenological profile will require 
Higgs studies at the LHC. 



15 



2.2.3 Standard Model Higgs Boson Searches at the LHC 



Production rates for the Higgs boson in the Standard Model are significantly larger at the LHC [ 44j,|69 |. 



The dominant Higgs production process, gluon fusion, can be exploited in conjunction with a variety of 
other channels, e.g., WW/ZZ fusion of the Higgs boson and Higgs radiation off top quarks [|7T| f7ll . f7ll . f^ . |rllJ7[)| 



Integrated luminosities between 30 and 100 fb , achievable within the first few years of LHC operation, 
will be sufficient to cover the entire canonical Higgs mass range of the Standard Model up to values 
close to 1 TeV with a significance greater than 5a. The required LHC luminosities for a Higgs discovery 
in various channels are shown in fig. 0(b) [72]. Thus, there is no escape route for the SM Higgs boson 
at the LHC. 

The properties of the SM Higgs boson can be determined with some precision at the LHC. The 
hsM ZZ^ — > £ + £~£ + £~ channel allows for an accurate Higgs mass determination of about 0.1% for 
120 GeV < rrih SM < 400 GeV, assuming an integrated luminosity of 300 fb" 1 |F(J. For larger Higgs 



masses, the precision in the Higgs mass measurement deteriorates due to the effect of the increasing Higgs 
width; nevertheless a 1% Higgs mass measurement is possible for m hsM ~ 700 GeV. The Higgs width 
can be extracted with a precision of 5 to 6% over the mass range 300 — 700 GeV from the Breit-Wigner 



shape of the Higgs resonance [76|]. Below 300 GeV, the instrumental resolution becomes larger than the 
Higgs width, and the accuracy of the Higgs width measurement degrades. For example, the four-lepton 
invariant mass spectrum from h$M — > ZZ yields a precision of about 25% at rrih SM = 240 GeV . For 
lower Higgs masses, indirect methods must be employed to measure the Higgs width. 

For Higgs masses below 200 GeV, a number of different Higgs decay channels can be studied at the 
LHC. The most relevant processes are 

99 -> h SM -> 77 , 
99 - h SM - VV& , 

qq -> qqV^V^ -> qqh SM , h SM -> 77, r + T~, VV& , 
gg, qq -> tih SM , h SM -> bb, 77, WW { * ] , 

where V = W or Z. The gluon-gluon fusion mechanism is the dominant Higgs production mechanism at 
the LHC, yielding a total cross-section of about 30 pb [15 pb] for rrih SM = 120 GeV [mh SM = 200 GeV]. 
One also has appreciable Higgs production via VV electroweak gauge boson fusion, with a total cross- 
section of about 6 pb [3 pb] for the Higgs masses quoted above. The electroweak gauge boson fusion 
mechanism can be separated from the gluon fusion process by employing a forward jet tag and central 
jet vetoing techniques. Note that for 2mw < m hsM < 2m z , the Higgs branching ratio to ZZ* is quite 
suppressed with respect to WW (since one of the Z bosons is off-shell). Hence, in this mass window, 
hsM — ► W + W~ — > £ + u£~u is the main Higgs discovery channel |65| , as exhibited in fig. 0(b). 

The cross-section for tih^M production can be significant for Higgs masses in the intermediate mass 



range [ iTf , 0.8 pb [0.2 pb] at rn,h SM = 120 GeV [mh SM = 200 GeV], although this cross-section falls faster 
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with Higgs mass as compared to the gluon and gauge boson fusion mechanisms. Finally, we note that 
the preferred channel at the Tevatron, qq — > Wh$M - > ivbb, is not a discovery mode at the LHC due 
in part to the larger background cross-sections at -^i = 14 TeV. Nevertheless, with 300 fb _1 of data, it 



may be possible to observe a SM Higgs signal (with S/yB > 5) if rrih SM < 125 GeV JfTj. 

The measurements of Higgs decay branching ratios can be used to infer the ratios of the Higgs 
couplings, and provide an important first step in clarifying the nature of the Higgs boson [78,79,30]. 
These can be extracted from a variety of Higgs signals that are observable over a limited range of Higgs 
masses. In the mass range 110 GeV < rrih SM < 150 GeV, the Higgs boson can be detected [with 100 fb _1 
of data] in the 77 and the r + r~ channels indicated above. (The fi + fi~ channel was considered in ref. [^0| . 
With 300 fb -1 of data, a 3a excess above background may be possible for 110 GeV < rrih SM < 140 GeV.) 
For m hsM > 130 GeV, the Higgs boson can also be detected in gluon-gluon fusion through its decay to 
WWW, with both final gauge bosons decaying leptonically |8T| , and to ZZ^*> in the four-lepton decay 
mode [|70| , |74| . There is additional sensitivity to Higgs production via VV fusion followed by its decay to 
WWW for rrih SM > 120 GeV. These data can be used to extract the ratios of the Higgs partial widths 
to gluon pairs, photon pairs, r + r~, and W + W~ |j78|JT9[] . The expected accuracies in Higgs width ratios, 
partial widths, and the total Higgs width are exhibited in fig. ||. These results are obtained under the 
assumption that the partial Higgs widths to W + W~ and ZZ are fixed by electroweak gauge invariance, 
and the ratio of the partial Higgs widths to bb and r + r~ are fixed by the universality of Higgs couplings 
to down-type fermions. One can then extract the total Higgs width under the assumption that all other 
unobserved modes, in the Standard Model and beyond, possess small branching ratios of order 1%. 
Finally, we note that the specific Lorentz structure predicted for the hg M W + W~ coupling by the Higgs 
mechanism can be tested in angular correlations between the spectator jets in WW fusion of the Higgs 
boson at the LHC [7[J. 

With an integrated luminosity of 100 fb _1 per experiment, the relative accuracy expected at the 
LHC for various ratios of Higgs partial widths Tj (or equivalently, ratios of Higgs squared-couplings) 
range from 10% to 30%, as shown in fig. |8|. The ratio T T /Tw measures the coupling of down-type 
fermions relative to the Higgs couplings to gauge bosons. To the extent that the one-loop /ism77 
amplitude is dominated by the W-loop, the partial width ratio r r /T 7 probes the same relationship. In 
contrast, under the usual assumption that the one- loop h^rngg amplitude is dominated by the top-quark 
loop, the ratio T g /Tw probes the coupling of up-type fermions relative to the HsmWW coupling. In 
addition, one can measure the hswiti coupling directly by making use of the tihsM production mode 
at the LHC. Recent studies suggest that for an integrated luminosity of 100 fb _1 , this signal is viable 



for the hsM —> bb and /i SM — > t + t~ [S3] decay modes if rrih SM ^ 130-140 GeV, and for the 

hsM — ¥ WW^ decay mode for 130 < mh SM < 200 GeV |55]. In this way, one expects to be able to 
measure the Higgs-top quark Yukawa coupling with a relative accuracy in the range of 10-20%. 

Finally, the measurement of the triple Higgs self-coupling with an accuracy of order 25% can be 
obtained, in a very limited Higgs mass window, with 3000 fb _1 of data [|86| (which requires at least a 
factor-of-ten luminosity upgrade of the LHC [87]). 
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Figure 8: Relative accuracy expected at the LHC with 200 fb _1 of data for (a) various ratios of Higgs boson partial widths 
and (b) the indirect determination of partial and total widths. Expectations for width ratios assume W, Z universality; 
indirect width measurements also assume b, r universality and a small branching ratio for unobserved modes. Taken from 
the parton-level analysis of ref. J78| . 



2.3 Standard Model Higgs Boson Searches at the LC 

The next generation of high energy e + e~ linear colliders is expected to operate at energies from 300 GeV 
up to about 1 TeV (JLC, NLC, TESLA), henceforth referred to as the LC []8~8,|89,90|. The possibility of a 
multi-TeV linear collider operating in an energy range of 3-5 TeV (CLIC) is also under study plj] . With 
the expected high luminosities up to 1 ab -1 , accumulated within a few years in a clean experimental 
environment, these colliders are ideal instruments for reconstructing the mechanism of electroweak 
symmetry breaking in a comprehensive and conclusive form. 

Weakly-coupled electroweak symmetry breaking dynamics involving an elementary scalar Higgs field 
can be established experimentally in three steps. First, the Higgs boson must be observed clearly and 
unambiguously, and its basic properties — mass, width, spin and C and P quantum numbers — must be 
determined. Second, the couplings of the Higgs boson to the and Z bosons and to leptons and 
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Figure 9: (a) The Higgs-strahlung and WW fusion production cross-sections as functions of rrih sm for s/s = 350 GeV, 



500 GeV and 800 GeV. (b)The cross-section for e + e~ — > tthsM, including NLO QCD corrections [ 102 1 , as a function of 
rrih SM for yfs = 500 GeV and 800 GeV with the expected experimental accuracy for m/i SM 
with error bar for an integrated luminosity of 1000 fb -1 . Taken from ref. p%. 
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quarks must be measured. Demonstrating that these couplings scale with the mass of the corresponding 
particle would provide critical support for the Higgs mechanism based on scalar dynamics as the agent 
responsible for generating the masses of the fundamental particles. Finally, the Higgs potential must 
be reconstructed by measuring the self-coupling of the Higgs field. The specific form of the potential 
shifts the ground state to a non-zero value, thereby providing the mechanism for electroweak symmetry 
breaking based on the self-interactions of scalar fields. Essential elements of this program can be realized 
at a high-luminosity e + e~ linear collider, and high-precision analyses of the Higgs boson are possible in 
these machines |^,|9|,|94|. 

The main production mechanisms of the SM Higgs boson at the LC are the Higgs-strahlung pro- 
p5|,p5p, e + e" — > Zh SM , and the WW fusion process |)6| e + e~ — > u e u e W*W* — > v e v e hsM- With an 



cess 



accumulated luminosity of 500 fb" 1 , about 10 5 Higgs bosons can be produced by Higgs-strahlung in the 
theoretically preferred intermediate mass range below 200 GeV. As y/s is increased, the cross-section 
for the Higgs-strahlung process decreases as s _1 and is dominant at low energies, while the cross-section 
for the WW fusion process grows as ln(s/m? ) and dominates at high energies (see ref. |J7| for a con- 



venient form for the corresponding cross-sections), as shown in fig. |9|(a).Q The ZZ fusion mechanism, 



8 The cross-sections shown in fig. |](a) are based on a tree-level computation. The complete one-loop electroweak 
corrections for e + e~ — > Zh^M have been computed 1^,0, whereas only the leading one- loop corrections to the vector 
boson fusion process are known |99| . 
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e + e~ — > e + e~Z*Z* —> e + e~hsM, also contributes to Higgs production, with a cross-section suppressed 
with respect to that of WW fusion by a factor of 16 cos 4 6V/[1 + (1 — 4:sm 2 9w) 2 } — 9.4. The cross- 
sections for the Higgs-strahlung and the WW fusion processes are shown in fig. |9|(a) for three values of 
s/s. The Higgs-strahlung process, e + e~ — > Zh$M, with Z — > £ + £~, offers a very distinctive signature. 
For y/s = 350 and 500 GeV and an integrated luminosity of 500 fb _1 , this ensures the observation of the 
SM Higgs boson up to the production kinematical limit independently of its decay ||92|| . At y/s = 500 
GeV, the Higgs-strahlung and the WW fusion processes have approximately the same cross-sections, 
0(50 fb), for 100 GeV < m hsM < 200 GeV. At ^ = 800 GeV with 500 fb" 1 of data, the analysis of 
ref. ||100|| suggests that a Higgs boson with mass up to about 650 GeV will be observable at the LC. 



Finally, the process e + e — > tthsM |pLQl|| yields a distinctive signature consisting of two W bosons and 



four 6-quark jets, and can be observed at the LC given sufficient energy and luminosity if the Higgs 
mass is not too large. The QCD-corrected cross-sections for this process [ |102| ] for y/s = 500 GeV and 
800 GeV are shown in fig. |9|(b). 

The phenomenological profile of the Higgs boson can be determined by precision measurements. 
For example, consider the case of frih SM = 120 GeV at the LC with y/s = 350 GeV and 500 fb" 1 of 
data |9~2]] . The spin and parity of the Higgs boson can be determined unambiguously from the steep 



onset of the excitation curve in Higgs-strahlung near the threshold (see fig. |10](a) ||103|| ) and the angular 



correlations in this process [ |104j| . By measuring final state angular distributions and various angular and 
polarization asymmetries, one can check whether the Higgs boson is a state of definite CP, or whether it 



exhibits CP- violating behavior in its production and/or decays | 105 |. The Higgs mass can be measured 
to an accuracy of 40 MeV by reconstructing the Higgs boson in Zh^M production and combining the 
results from the various final state channels. The Higgs width can be inferred in a model-independent 
way, with an accuracy of about 6%, by combining the partial width to W + W~, accessible in the vector 
boson fusion process, with the WW* decay branching ratio. Similar results (with precisions within 
a factor of two of those quoted above) are obtained for larger Higgs masses in the intermediate mass 
regime. The hsu_ZZ, and hsMtt couplings can be determined (with some sensitivity to possible 

anomalous couplings, if present) by the optimal-observable method ||106| , which makes optimal use of 
the polarized angular distributions and asymmetries for Higgs-strahlung and hsMti production. 



Higgs decay branching ratios can be measured very precisely for rrih SM ^ 150 GeV [|107| , |108| , |109| . |110 |. 
When such measurements are combined with measurements of Higgs production cross-sections, the 
absolute values of the Higgs couplings to the W ± and Z gauge bosons and the Yukawa couplings to 
leptons and quarks can be determined to a few percent in a model-independent way. In addition, 
the Higgs-top quark Yukawa coupling can be inferred from the cross-section for Higgs emission off 
ti pairs |102| , |111| . As an example, Table [I] exhibits the anticipated fractional uncertainties in the 
measurements of Higgs branching ratios for m>h SM = 120 GeV at the LC.[] Using this data, a program 



9 Here, BR(/igM - * gg) is assumed to be roughly equal to the Higgs branching ratio into light hadrons (i.e., excluding 
hadrons that contain c and b quarks). 
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Figure 10: (a) Simulated measurement of the e + e~ — > ZhsM cross-section for irih SM = 120 GeV with 20 fb _1 per 
point at three center of mass energies compared to the predictions for spin-0 (solid line) and typical examples of spin-1 
(dashed line) and spin-2 (dotted line) |105( |. (b) Cross-section for the double Higgs-strahlung process e + e~ — > ZhsuhsM 
at yfs = 500 GeV (solid line) and 800 GeV (dashed line) Jll4| . The data points show the accuracy for 1 ab _1 . 
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HFITTER was developed in ref. [107] to perform a Standard Model global fit based on the measurements 
of the ZhsM, vvhsM and tthsu cross-sections and the Higgs branching ratios listed in Table |l[ The 
output of the program is a set of Higgs couplings along with their fractional uncertainties (which are 
also exhibited in Table [I]). These results should be considered representative of what can eventually 
be achieved at the LC. For example, a comprehensive analysis of LC Higgs data will have to take 
into account the relevant electroweak loop corrections (which are not presently included in HFITTER). 
One should also note that theoretical uncertainties for the predicted Higgs couplings have not been 
taken into account in this analysis. The theoretical uncertainty in gh SM cc is the most significant among 
the channels listed in Table [l], due to the uncertainties in c quark-mass and in a s (which governs the 
running of the quark masses from the quark mass to the Higgs mass). Ref. ||112|| estimates a theoretical 
fractional uncertainty in g^ cZ of about 12%, significantly greater than the experimental uncertainty listed 
in Table [I]. In contrast, the theoretical fractional uncertainty in g hb i is about 1.8% due to the uncertainty 
of the 6-quark mass. Although this is less than the anticipated experimental uncertainty, it should not 
be neglected in the determination of the overall gh SM bb uncertainty. The theoretical uncertainties in the 
other channels listed above are not significant compared to the quoted experimental uncertainties. 

The measurement of the Higgs self-couplings is a very ambitious task that requires the highest 
luminosities possible at the LC, which possess unique capabilities for addressing this question. The 
trilinear Higgs self-coupling can be measured in double Higgs-strahlung, in which a virtual Higgs boson 
splits into two real Higgs particles in the final state [ 113 |. A simulation based on 1 ab^ 1 of data is 



exhibited in fig. |lO|(b) [ |114j| . In this way, for rrih SM = 120 GeV, the cubic term of the scalar potential can 
be established at the LC with a precision of about 20% [ |114| , |115| |. Such a measurement is a prerequisite 
for determining the form of the Higgs potential that is responsible for spontaneous electroweak symmetry 
breaking generated by scalar sector dynamics. Finally, the total SM Higgs width can be obtained 
indirectly by using r to t = ^h SM ww /BR(/ism — > WW*). The partial width is proportional to gl SMWW , 
so the fractional uncertainty in r tot can be obtained from the results of Table |l|. For mh SM < 150 GeV, 
an accuracy in the range of 5-10% can be achieved for the total Higgs width [ |116|| . 

If the SM Higgs mass is above 150 GeV, then the precision determination of Higgs couplings will 
have to be reconsidered. The Higgs branching ratios into cc, gg and r + r _ are now too small to be 
accurately measured. Due to the growing importance of the WW and ZZ modes, one can perform 
a precision measurement of the Higgs branching ratios to WW and ZZ, while the precision of the bb 
branching ratio is significantly reduced as rrih SM increases to 200 GeV and beyond ||117]| . Moreover, due 
to the rapid decline of tihsM production cross-section with increasing rrih SM , the Higgs-top quark Yukawa 
coupling cannot be extracted until mh SM > 2m t , at which point the tt/isM coupling can be obtained by 
observing Higgs bosons produced by vector boson fusion which subsequently decay to ti. The analysis of 
Ref. ||118|| finds that at the LC with a/s = 800 GeV and 1 ab" 1 of data, the tihsM Yukawa coupling can 



be determined with an accuracy of about 10% for a Higgs mass in the range 350 — 500 GeV. The total 
Higgs width can be obtained directly from measuring the Higgs boson line-shape if m hsM > 200 GeV. 
The e + e~ linear collider with center-of-mass energy ^/s can also be designed to operate in a 77 col- 
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lision mode. This is achieved by using Compton backscattered photons in the scattering of intense laser 
photons on the initial polarized e* 1 beams ||119| , |120|| . The resulting 77 center of mass energy is peaked 
for proper choices of machine parameters at about 0.8-y/i. The luminosity achievable as a function of the 
photon beam energy depends strongly on the machine parameters (in particular, the choice of laser polar- 
izations). The 77 collider provides additional opportunities for Higgs physics [ |1 20|J12 1|J1 22| , |1 23| , |1 24|| . The 
Higgs boson can be produced as an s-channel resonance in 77 collisions, and one can perform indepen- 
dent measurements of various Higgs couplings. For example, the product T(/ism — ► 77)BR(/ism —> bb) 
can be measured with a statistical accuracy of about 2 — 10% for 120 GeV< 77i^ SM < 160 GeV with about 
50 fb" 1 of data [ 122 , 123]J124j| . In order to reach such a precision, it is critical to control the overwhelm- 
ing two-jet background (with efficient 6-tagging) and overcome the irreducible 77 — > bb background by 
optimal use of the polarization of the photon beams and by judicious kinematic cuts. Knowledge of the 
QCD corrections to signal and background processes is essential for this task ||122 , [T25|| . 

Using values for BR(/igM — > bb) and BR(/igM — > 77) measured at the e + e~ linear collider, one can 
obtain a value for the total Higgs width with an error dominated by the expected error in BR(/ism — > 77). 
For heavier Higgs bosons, mh SM > 200 GeV, the total Higgs width can in principle be measured directly 
by tuning the collider to scan across the Higgs resonance. One can also use the polarization of the 
photon beams to measure various asymmetries in Higgs production and decay, which are sensitive to 



the CP quantum numbers of the Higgs boson [123 



3 The Higgs Bosons of Low-Energy Super symmetry 

Electroweak symmetry breaking dynamics driven by a weakly-coupled elementary scalar sector requires 
a mechanism for the stability of the electroweak symmetry breaking scale with respect to the Planck 
scale ||. Supersymmetry-breaking effects, whose origins may lie at energy scales much larger than 
1 TeV, can induce a radiative breaking of the electroweak symmetry due to the effects of the large 



Higgs-top quark Yukawa coupling | 126 |. In this way, the origin of the electroweak symmetry breaking 
scale is intimately tied to the mechanism of supersymmetry breaking. Thus, supersymmetry provides 
an explanation for the stability of the hierarchy of scales, provided that supersymmetry-breaking masses 
in the low-energy effective electroweak theory are of (9(1 TeV) or less ||. 

A fundamental theory of supersymmetry-breaking is presently unknown. Nevertheless, one can 
parameterize the low-energy theory in terms of the most general set of soft-supersymmetry-breaking 



terms ||127|| . The simplest realistic model of low-energy supersymmetry is a minimal supersymmetric 
extension of the Standard Model (MSSM), which employs the minimal supersymmetric particle spec- 
trum. However, even in this minimal model with the most general set of soft-supersymmetry-breaking 
terms, more than 100 new supersymmetric parameters are introduced [|1 28|| . Fortunately, most of these 
parameters have no impact on Higgs phenomenology. Thus, we will focus primarily on the Higgs sector 
of the MSSM and identify the parameters that govern the main properties of the Higgs bosons. 
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3.1 The Tree-Level Higgs Sector of the MSSM 

Both hypercharge Y — — 1 and Y = +1 complex Higgs doublets are required in any Higgs sector 
of an anomaly-free supersymmetric extension of the Standard Model. The supersymmetric structure 
of the theory also requires (at least) two Higgs doublets to generate mass for both "up"-type and 
"down" -type quarks (and charged leptons) ||. Thus, the MSSM contains the particle spectrum of a 
two-Higgs-doublet extension of the Standard Model and the corresponding supersymmetric partners. 

The two-doublet Higgs sector ||129|| contains eight scalar degrees of freedom: one complex Y = — 1 
doublet, < 3?d= (^d^d) an d one complex Y = +1 doublet, <& u = (*£« The notation reflects the 

form of the MSSM Higgs sector coupling to fermions: $° [$°] couples exclusively to down-type [up- 
type] fermion pairs. When the Higgs potential is minimized, the neutral components of the Higgs fields 
acquire vacuum expectation valuesQ 



($ d) =4= , <*„>=4= : , (5) 





where the normalization has been chosen such that v 2 = v^ + v 2 = ^m^j g 2 = (246 GeV) 2 . Spontaneous 
electroweak symmetry breaking results in three Goldstone bosons, which are absorbed and become the 
longitudinal components of the and Z. The remaining five physical Higgs particles consist of a 
charged Higgs pair 

H ± = $ffn.n/3 + $±coa/3, (6) 



one CP-odd scalar 

and two CP-even scalars: 



A = V2 flm $° sin /? + Im $° cos (3j , (7) 



h = - ( V2 Re $° - v d ) sin a + ( V2 Re $° - v u ) cos a , 

H = (v / 2Re$S-^)cosa + (v / 2Re$°-i; u )sina, (8) 

(with rrih < mu)- The angle a arises when the CP-even Higgs squared-mass matrix (in the $° — $° 
basis) is diagonalized to obtain the physical CP-even Higgs states (explicit formulae will be given below). 

The supersymmetric structure of the theory imposes constraints on the Higgs sector of the model. 
For example, the Higgs self-interactions are not independent parameters; they can be expressed in terms 
of the electroweak gauge coupling constants. As a result, all Higgs sector parameters at tree-level are 
determined by two free parameters: the ratio of the two neutral Higgs field vacuum expectation values, 

tan/5 = —, (9) 

v d 

10 The phases of the Higgs fields can be chosen such that the vacuum expectation values are real and positive. That 
is, the tree-level MSSM Higgs sector conserves CP, which implies that the neutral Higgs mass eigenstates possess definite 
CP quantum numbers. 
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and one Higgs mass, conveniently chosen to be tua- In particular, 

m 2 H± =m 2 A + m 2 w , (10) 
and the CP-even Higgs bosons h and H are eigenstates of the following squared-mass matrix 



m\ sin 2 /3 + m% cos 2 (3 -{m\ + m|) sin (3 cos (3 , 
V — { m \ + m z) sin/3 cos/3 m 2 ^ cos 2 /3 + m| sin 2 /3 

The eigenvalues of M.1 are the squared- masses of the two CP-even Higgs scalars 

m H,h = \ (™ 2 a + m% ± \] [m\ + m 2 z ) 2 - Am 2 z m 2 A cos 2 2/3 ) , (12) 

and a is the angle that diagonalizes the CP-even Higgs squared-mass matrix. From the above results, 
one obtains: 

cos 2 ((3-a)= m l {m l~ ml l . (13) 
m A ( m H ~ m h) 

In the convention where tan/3 is positive (i.e., < (3 < vr/2), the angle a lies in the range —tt/2 < a < 0. 

An important consequence of eq. fll2|) is that there is an upper bound to the mass of the light 
CP-even Higgs boson, h. One finds that: 

m h < m z| cos2/3| < m z ■ (14) 

This is in marked contrast to the Standard Model, in which the theory does not constrain the value of 
m/i SM at tree-level. The origin of this difference is easy to ascertain. In the Standard Model, Tn^ = |At> 2 
is proportional to the Higgs self-coupling A, which is a free parameter. On the other hand, all Higgs 
self-coupling parameters of the MSSM are related to the squares of the electroweak gauge couplings. 

Note that the Higgs mass inequality [eq. (0)] is saturated in the limit of large rriA- In the limit of 
vn a ^> friz, the expressions for the Higgs masses and mixing angle simplify and one finds 

m 2 h ~ m| cos 2 2(3, (15) 



x\ ~ m 2 A + m\ - 

,2 _2 i _2 



m\ + m 2 z sin 2 2(3 , (16) 



m H ± = m A + m w , (17) 

9 . „ . mi sin 2 4/3 . , 

cos 2 (/3-a)^ Z 4m 4 ~ • (18) 

Two consequences are immediately apparent. First, rriA — tuh — vnu±, up to corrections of 0(m 2 z /mA)- 
Second, cos(/3 — a) = up to corrections of 0(m 2 z /m 2 A ). This limit is known as the decoupling limit [ 180 
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because when is large, there exists an effective low-energy theory below the scale of in which the 
effective Higgs sector consists only of one CP-even Higgs boson, h. As we shall demonstrate below, the 
tree-level couplings of h are precisely those of the Standard Model Higgs boson when cos(/? — a) = 0. 
From eq. (|T8|), one can also derive: 

cota = -tan/5 f tan p cos 2/3 + O — f . (19) 

This result will prove useful in evaluating the CP-even Higgs boson couplings to fermion pairs in the 
decoupling limit. 

The phenomenology of the Higgs sector depends in detail on the various couplings of the Higgs 
bosons to gauge bosons, Higgs bosons and fermions. The couplings of the two CP-even Higgs bosons 
to W and Z pairs are given in terms of the angles a and f3 by 

9hW = 9v m v sin (^ - a ) > 9 H vv = 9v m v cos (P ~ a ) > ( 20 ) 

where g v = Imyjv for V = W or Z. There are no tree-level couplings of A or H ± to VV. The 
couplings of V to two neutral Higgs bosons (which must have opposite CP-quantum numbers) are given 
by g<i>Az(p<t> — Pa), where (ft = h or H and the momenta and pa point into the vertex, and 

gcos((3- a) -gsm((3 - a) 

9hAZ ~ 2 C os9 w ' 9haz ~ 2cdb9 w ' 1 j 

From the expressions above, we see that the following sum rules must hold separately for V = W and Z: 

9hw + 9lw = 9v m2 v > ( 22 ) 

2 2 9 2 

9hAz + 9haz = A „„„ 9 /i ' 

(23) 



4 COS 2 
2 2 

4^ + 4m|^ = -^-, = /i ?J ff. (24) 
cos z v w 



Similar considerations also hold for the coupling of h and H to W^H^. Four-point couplings of vector 
bosons and Higgs bosons can be found in ref. 0. The properties of the three-point and four-point Higgs 
boson- vector boson couplings are conveniently summarized by listing the couplings that are proportional 
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to either sin(/3 — a) or cos(/3 — a), and the couplings that are independent of a and (3 0: 



cos(/3 — a) 

HW + W~ 

HZZ 

ZAh 

W ± H T h 

ZW ± H T h 

1 W ± H^h 



sin(/3 — at) 

hW+W~ 

hZZ 

ZAH 

W ± H^H 

ZW ± H T H 



angle- independent 



ZH + H~, -fH + H~ 
W ± H^A 
ZW ± H^A 
1 W ± H^A 

VV(j)(j>, VVAA, VVH + H~ 



(25) 



where <p = h or H and VV = W + W~ , ZZ, Zj or 77. Note in particular that all vertices in the theory 
that contain at least one vector boson and exactly one non-minimal Higgs boson state (H, A or H^) 
are proportional to cos(/3 — a). This can be understood as a consequence of unitarity sum rules which 
must be satisfied by the tree-level amplitudes of the theory [ ^j25lJT3l| . |132 l. 

In the MSSM, the tree-level Higgs couplings to fermions obey the following property: $° couples 
exclusively to down-type fermion pairs and $° couples exclusively to up-type fermion pairs. This pattern 
of Higgs-fermion couplings defines the Type-II two-Higgs-doublet model ||133| j2]1. The gauge- invariant 
Type-II Yukawa interactions (using 3rd family notation) are given by: 



C 



Yukawa 



tP L t$° u - iP L b<5>+] + h b \bP L b$° d - bP L t$a] + h.c. 



(26) 



where Pl 



|(1 — 75) is the left-handed projection operator. [Note that (^ 1 P L ^/ 2 y = where 
- R = -yi -r 75).] Fermion masses are generated when the neutral Higgs components acquire vacuum 
expectation values. Inserting eq. (§) into eq. ( ^6|) yields a relation between the quark masses and the 
Yukawa couplings: 

V2m b y/2m b V2m t V2m t 
h = = > h t = = — ^7 . (27) 



v cos (3 



v sin (3 



Similarly, one can define the Yukawa coupling of the Higgs boson to r-leptons (the latter is a down-type 
fermion). The couplings of the physical Higgs bosons to the third generation fermions is obtained from 
eq. fl2HD by using eqs. 111 particular, the couplings of the neutral Higgs bosons to // relative to 

the Standard Model value, gmf/2mw, are given by 



hbb (or hr + T~ 



htt 



sin a 

cos (3 
cos a 
sin (3 



= sin(/3 — at) — tan (3 cos(/3 — a) 
sin(/5 — a) + cot (3 cos((3 — a) , 



(28) 
(29) 



27 



— cos en 

Hbb (or Ht + t~): = cos(/3 - a) + tan (3 sin(/3 - a) , (30) 

cos p 

sin a 

Htt : = cos(/3 — a) — cot (3 sin(/3 — a) , (31) 

sin p 

Abb (or Ar + r") : 75 tan /3, (32) 

Att : 75 cot /3 , (33) 

(the 75 indicates a pseudoscalar coupling), and the charged Higgs boson couplings to fermion pairs, 
with all particles pointing into the vertex, are given by 



9H~tb 



y/2m w 
9h~ + = 9 



m t cot (3 Pr + nib tan (3 Pl 



m T tan f3 Pl 



(34) 



(35) 



We next examine the behavior of the Higgs couplings at large tan/3. This limit is of particular 
interest since at large tan (3, some of the Higgs couplings to down-type fermions can be significantly 
enhanced.^ Consider two large tan (3 regions of interest: (i) If ijia ^> mz, then the decoupling limit is 
reached, in which | cos(/3 — a)| <C 1 and mn — ttia- From eqs. (|T5[) — (P3[) , it follows that the bbH and bbA 
couplings have equal strength and are significantly enhanced (by a factor of tan/3) relative to the bbh^m 
coupling, whereas the VVH coupling is negligibly small. In contrast, the values of the VVh and bbh 
couplings are equal to the corresponding couplings of the Standard Model Higgs boson. To show that 



the value of the bbh coupling [eq. (£8|)] reduces to that of bbh^M in the decoupling limit, note that eq. ( [18]) 
implies that | tan f3 cos(/3 — a)\ 1 when wla 3> mz even when tan (3 ^> 1. Indeed, h is a SM-like Higgs 
boson, (ii) If rriA ^ rn z and tan/3 ^> 1, then | cos(/3 — a)\ ~ 1 [see fig. [11]] and m h ~ m^. In this case, 
the bbh and bbA couplings have equal strength and are significantly enhanced (by a factor of tan /3) 
relative to the bbhsM coupling, while the VVh coupling is negligibly small. Using eq. (|22"D it follows that 
the VVH coupling is equal in strength to the VVhsM coupling. In this case, it is conventional to refer 
to if as a SM-like Higgs boson. However, this nomenclature is somewhat inaccurate, since the value of 
the bbH coupling can differ from the corresponding bbhsM coupling when tan/3 3> 1 [since in case (ii), 
where | sin(/3 — a)\ 1, the product tan f3 sin(/3 — a) need not be particularly small]. Note that in both 
cases (i) and ii) above, only two of the three neutral Higgs bosons have enhanced couplings to bb. 

The decoupling limit of tua 3> rn>z is effective for all values of tan/3. It is easy to check that the 
pattern of all Higgs couplings displayed in eqs. (f20[)-(|33D respect the decoupling limit. That is, in the 



models of low-energy supersymmetry, there is some theoretical prejudice that suggests that 1 < tan/3 < mt/rrib, 
with the fermion running masses evaluated at the electroweak scale. For example, tan (3 < 1 [tan (3 > m t / rrib] is disfavored 
since in this case, the Higgs-top-quark [Higgs-bottom-quark] Yukawa coupling blows up at an energy scale significantly 
below the Planck scale. 
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Figure 11: The value of cos 2 (/3 — a) is shown as a function of tua for two choices of tan/3 = 3 and tan/3 = 30. When 
radiative-corrections are included, one can define an approximate loop-corrected angle a as a function of ttia, tan/3 and 
the MSSM parameters. In the figures above, we have incorporated radiative corrections, assuming that Msusy = Mq = 
Mjj = Mu = 1 TcV. In addition, two extreme cases for the squark mixing parameters arc shown (see Sections [T^ and 
3.3 for further discussion of the radiative corrections and their dependence on the supersymmetric parameters). The 
decoupling effect expected from eq. ([T8|), in which cos 2 (/3 — a) cx m\jm A for itia 3> rnz, continues to hold even when 
radiative corrections are included. 



limit where 3> rn>z, cos(/3 — a) = 0{m 2 z /m 2 A ), which means that the h couplings to Standard Model 
particles approach values corresponding precisely to the couplings of the SM Higgs boson. There is a 
significant region of MSSM Higgs sector parameter space in which the decoupling limit applies, because 



cos(/3 — a) approaches zero quite rapidly once is larger than about 200 GeV, as shown in fig. |TT 



As a result, over a significant region of the MSSM parameter space, the search for the lightest CP-even 
Higgs boson of the MSSM is equivalent to the search for the Standard Model Higgs boson. This result 
is more general; in many theories of non-minimal Higgs sectors, there is a significant portion of the 
parameter space that approximates the decoupling limit. Consequently, simulations of the Standard 
Model Higgs signal are also relevant for exploring the more general Higgs sector. 



3.2 Radiatively-Corrected MSSM Higgs Masses 



The discussion of Section |3.1| was based on a tree-level analysis of the Higgs sector. However, radiative 
corrections can have a significant impact on the predicted values of Higgs masses and couplings. The 
radiative corrections involve both loops of Standard Model particles and loops of supersymmetric part- 
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ners. The dominant effects arise from loops involving the third generation quarks and squarks and are 
proportional to the corresponding Yukawa couplings. Thus, we first review the parameters that con- 
trol the masses and mixing of the third-generation squarks. (We shall neglect intergenerational mixing 
effects, which have little impact on the discussion that follows.) 

For each left-handed and right-handed quark of fixed flavor, q, there is a corresponding supersymmet- 
ric partner and cfn, respectively. These are the so-called interaction eigenstates, which mix according 
to the squark squared-mass matrix. The mixing angle that diagonalizes the squark mass matrix will be 
denoted by 9~. The squark mass eigenstates, denoted by q\ and q 2 , are obtained by diagonalizing the 
following 2x2 matrix 

Ml + m 2 f + D L m f X f \ 

, (36) 

mfXf M 2 R + m) + D R J 

where Dj = (T 3 f — e/ sin 2 6w)m 2 z cos 2/? and Dr = e/sin 2 9wm 2 z cos 2/3. In addition, f = t, Mr = Mjj, 
e t = 2/3 and T 3 f = 1/2 for the top-squark squared-mass matrix, and / = b, Mr = Mp, e& = —1/3 and 
T 2 f = —1/2 for the bottom-squark squared-mass matrix. The squark mixing parameters are given by 

X t = A t — n cot (3 , X b = A b — ji tan (3 . (37) 

Thus, the top-squark and bottom-squark masses and mixing angles depend on the supersymmetric 
Higgsino mass parameter \x and the soft-supersymmetry-breaking parameters: Mq, Mjj, Mjj, A t and 
Ab. For simplicity, we shall initially assume that A t , A b and \i are real parameters. That is, we neglect 
possible CP-violating effects that can enter the MSSM Higgs sector via radiative corrections. The 



impact on new MSSM sources of CP-violation on the Higgs sector will be addressed in Section p. 2. 3 



3.2.1 Radiatively-corrected Higgs masses in the CP-conserving MSSM 

The radiative corrections to the Higgs squared-masses have been computed by a number of techniques, 
and using a variety of approximations such as the effective potential at one-loop ||13 1| . |135| . |136| and 



two-loops |T37| , pg|CT , [Hqi , and diagrammatic methods PTjP^|H3l|T~ni . |l ^[gg . Complete one- 



loop diagrammatic computations of the MSSM Higgs masses have been presented by a number of 
groups |p.43| , |144| ; and partial two- loop diagrammatic results are also known ||145| , |146|| . These in- 



clude the 0(m 2 h 2 a s ) contributions to the neutral CP-even Higgs boson squared-masses in the on-shell 
scheme | |146| |. Finally, renormalization group methods (to be discussed further below) provide a pow- 



erful technique for identifying many of the most important contributions to the radiatively corrected 
Higgs masses [ |147| . |148| . |149|jl50[ | . Typical results for the radiatively corrected value of rrih as a function 
of the relevant supersymmetric parameters are shown in fig. [1^. 

One of the most striking effects of the radiative corrections to the MSSM Higgs sector is the mod- 



ification of the upper bound of the light CP-even Higgs mass, as first noted in refs. ||134|| and ||141 



Consider the region of parameter space where tan/5 is large and 3> rriz- In this limit, the tree-level 
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Figure 1 2 : The radiatively corrected light CP-even Higgs mass is plotted (a) as a function of X t , where X t = At — n cot /3, 
for M t = 174.3 GeV and two choices of tan/3 = 3 and 30, and (b) as a function of tan/3, for the maximal mixing [upper 
band] and minimal mixing [lower band] benchmark cases. In (b), the central value of the shaded bands corresponds to 
Mt = 175 GeV, while the upper [lower] edge of the bands correspond to increasing [decreasing] Mt by 5 GeV. In both 
(a) and (b), via = 1 TeV and the diagonal soft squark squared-masses are assumed to be degenerate: Msusy = Mq = 
M v = M D = 1 TeV. 



prediction for corresponds to its theoretical upper bound, = rriz- Including radiative corrections, 
the theoretical upper bound is increased. The dominant effect arises from an incomplete cancellation^ 
of the top-quark and top-squark loops (these effects cancel in the exact supersymmetric limit). The 
qualitative behavior of the radiative corrections can be most easily seen in the large top squark mass 
limit, where in addition, the splitting of the two diagonal entries and the off-diagonal entry of the 
top-squark squared-mass matrix are both small in comparison to the average of the two top-squark 
squared-masses: 

i(Mg+M|). (38) 



Ml 



In this case, the upper bound on the lightest CP-even Higgs mass is approximately given by 



m 



< 



rrir 



3g 2 mj 




12M 2 S/ 



(39) 



12 In certain regions of parameter space (corresponding to large tan /3 and large values of /i), the incomplete cancellation 
of the bottom-quark and bottom-squark loops can be as important as the corresponding top sector contributions. For 
simplicity, we ignore this contribution in eq. (fey). 
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The more complete treatments of the radiative corrections cited above show that eq. (^) somewhat 
overestimates the true upper bound of rrih- Nevertheless, eq. ( p9[) correctly reflects some noteworthy 
features of the more precise result. First, the increase of the light CP-even Higgs mass bound beyond 
rriz can be significant. This is a consequence of the mf enhancement of the one- loop radiative correction. 
Second, the dependence of the light Higgs mass on the top-squark mixing parameter X t implies that (for 
a given value of Ms) the upper bound of the light Higgs mass initially increases with X t and reaches 
its maximal value for X t = y/EMg- This point is referred to as the maximal mixing case (whereas 
Xt = corresponds to the minimal mixing case). In a more complete computation that includes both 
two-loop logarithmic and non-logarithmic corrections, the X t values corresponding to maximal and 
minimal mixing are shifted and exhibit an asymmetry under X t — » —X t as shown in fig. In the 
numerical analysis presented in this and subsequent figures in this section, we assume for simplicity 
that the third generation diagonal soft-supersymmetry-breaking squark squared-masses are degenerate: 
M S usy = Mq = Mjj = M D , which defines the parameter M SU sy-0 

Third, note the logarithmic sensitivity to the top-squark masses. Naturalness arguments that un- 
derlie low-energy supersymmetry imply that the supersymmetric particle masses should not be larger 
than a few TeV. Still, the precise upper bound on the light Higgs mass depends on the specific choice 
for the upper limit of the top-squark masses. The dependence of the light Higgs mass obtained by the 
more complete computation as a function of Msusy is shown in fig. [T^Q 



As noted above, the largest contribution to the one-loop radiative corrections is enhanced by a factor 
of m\ and grows logarithmically with the top squark mass. Thus, higher order radiative corrections can 
be non-negligible for large top squark masses, in which case the large logarithms must be resummed. 
Renormalization group (RG) techniques for resumming the leading logarithms have been developed by 
a number of authors | |147| , |14£| , |i4iJ[ | . The computation of the RG-improved one- loop corrections requires 
numerical integration of a coupled set of RG equations ||148|| . Although this procedure has been carried 



out in the literature, the analysis is unwieldy and not easily amenable to large-scale Monte-Carlo studies. 
It turns out that over most of the parameter range, it is sufficient to include the leading and sub-leading 
logarithms at two-loop order. (Some additional non- logarithmic terms, which cannot be ascertained 
by the renormalization group method, must also be included ||151|| .) Compact analytic expressions 
have been obtained for the dominant one and two-loop contributions to the matrix elements of the 
radiatively-corrected CP-even Higgs squared-mass matrix: 

/Ml M\ 2 \ 

M 2 = I ^ l 2 )=Ml + 5M\ (40) 
\M~ 12 M 2 22 ) 

where the tree-level contribution Ai^ was given in eq. (]TT|) and 5Ai 2 is the contribution from the 



13 We also assume that AigusY 3> mt, in which case it follows that Mj ~ M| USY up to corrections of C(mj /Af| ugY ). 
14 The flattening of the curves in fig. [Ds] as a function of Msusy in the maximal mixing scenario is due to the squark- 
mixing contributions at two- loops which partially cancel the contributions that grow logarithmically with Msusy- 
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Figure 13: The radiatively corrected light CP-even Higgs mass is plotted as a function of Msusy = Mq = Mu = Md, 
for Mt = 174.3 GeV, uia = 1 TeV and two choices of tan/3 = 3 and tan/? = 30. Maximal mixing and minimal mixing 
are defined according to the value of X t that yields the maximal and minimal Higgs mass as shown in fig. |l^(a). 



radiative corrections. The dominant corrections to Ai 2 , coming from the one-loop top and bottom 
quark and top and bottom squark contributions plus the two-loop leading logarithmic contributions, 
are given to 0{h\, h\) by [p.49|J15q , [T52 
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where sp = sin/?, cp = cos/5, and the coefficients Qj are: 



Cij — 



32tt 2 
33 



(44) 



(t n ,t 12 ,t 21 ,t 22 ,t 31 ,t 32 ) = (12,-4,6,-10,9,-7) and (b n ,b 12 ,b 21 ,b 22 ,b 31 ,b 32 ) = (-4,12,2,18,-1,15). 
Above, h t and h b are the top and bottom quark Yukawa couplings [see eqs. (p5|)-([57|)1, g 3 is the strong 
QCD coupling, v = 246 GeV is the SM Higgs vacuum expectation value, and Mj = ^(M? + M~) is 
the average squared top squark mass.0 The SAifj also depend on the MSSM parameters A t , A b and 
ji that enter the off-diagonal top-squark and bottom-squark squared-mass matrices. We employ the 
following notation: p, = fi/M s , a t = A t /M s , a b = A b /M s and x t = X t /M s , where X t = A t — /icot/3. 
Diagonalizing the CP-even Higgs squared-mass matrix yields radiatively-corrected values for m\, m\ 
and the mixing angle The end result is a prediction for the Higgs mass in terms of running 

parameters in the MS scheme. It is a simple matter to relate these parameters to the corresponding 
on-shell parameters used in the diagrammatic calculations ||138| , |151|| . 

Additional non-logarithmic two-loop contributions, which can generate a non-negligible shift in the 
Higgs mass (of a few GeV), must also be included.^ A compact analytical expression that incorporates 
these effects at 0(mfhfa s ) was given in ref. ||154|| (with further refinements provided by ref. ||155|| to take 
into account the possibility of arbitrary top-squark splitting), and the corresponding corrections pro- 
portional to hla s can be found in ref. ||140|| . An important source of such contributions are the one- loop 
super symmetric threshold corrections to the relation between the Higgs-top-quark and Higgs-bottom- 
quark Yukawa couplings and the corresponding quark masses [eqs. fl56|) and fl57D]. These generate a 
non-logarithmic two-loop shift of the radiatively corrected Higgs mass proportional to the correspond- 
ing squark mixing parameters. One consequence of these contributions [ |151|| is the asymmetry in the 
predicted value of m h under X t — > —X t as noted in fig. 0(a). Recently, the computation of m h has 
been further refined by the inclusion of genuine two-loop corrections of 0(m^hj) ||139|| , and estimates 



of the two-loop corrections proportional to hf/jf and h\ [ |140|| (which can be numerically relevant for 
values of tan/3 > m t /m b ). These non-logarithmic corrections, which depend on the third generation 
squark mixing parameters, can slightly increase the value of the radiatively-corrected Higgs mass. 

The numerical results displayed in figs. |TTHT4l are based on the calculations of refs. ||149|| and ||1 50 . 
with improvements as described in refs. ||146|| and [|151[] . The supersymmetric parameters in the maximal 
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and minimal mixing cases have been chosen according to the first two benchmark scenarios of ref. 
Of particular interest is the upper bound for the lightest CP-even Higgs mass (mh). At fixed tan/3, the 
maximal value of m,h is reached for vtla ^> raz (see fig. |14]). Taking itla large, fig. [L2|(b) illustrates that 
the maximal value of the lightest CP-even Higgs mass bound is realized at large tan/3 in the case of 



15 Eqs. ([tl|)-([t3|) have been derived under the assumption that |M? 



-M?\/(Ml 



' ) <C 1. The approximate forms of 

eqs. fl4l|) -([43|) are sufficient to provide insight on the dependence of the radiatively-corrected Higgs masses and couplings 
on the MSSM parameters, although our numerical work is based on more exact forms for these expressions. 

16 Although M.\ 2 is negative at tree level (implying that — 7r/2 < a < 0), it is possible that radiative corrections flip the 

-tt/2 < a < tt/2. 



sign of M.\ 2 - Thus, the range of the radiatively corrected angle a can be taken to be 



17 An improved procedure for computing the radiatively-corrected neutral Higgs mass matrix and the charged Higgs mass 
in a self-consistent way (including possible CP-violating effects), which incorporates one- loop supersymmetric threshold 



corrections to the Higgs-top-quark and Higgs-bottom-quark Yukawa couplings, can be found in ref. |153| 



34 



> 

1) 



(0 



250 



200 



150 



100 




_L 



tan = 3 

_^ , L_ 



100 



150 
m 4 



200 



250 



(GeV) 



Figure 14: Lightest CP-even Higgs mass (m^), heaviest CP-even Higgs mass (ma) and charged Higgs mass (m^-±) as a 
function of mj for two choices of tan/3 = 3 and tan/3 = 30. Here, we have taken M t = 174.3 GeV, and we have assumed 
that the diagonal soft squark squared-masses are degenerate: Msusy = Mq = Mu = Mjj = 1 TcV. In addition, we 
choose the other supcrsymmetric parameters corresponding to the maximal mixing scenario. The slight increase in the 
charged Higgs mass as tan/3 is increased from 3 to 30 is a consequence of the radiative corrections. 



maximal mixing. For each value of tan/3, we denote the maximum value of mh by m™ ax (tan/5) [this 
value also depends on the third-generation squark mixing parameters]. Allowing for the uncertainty 
in the measured value of m t and the uncertainty inherent in the theoretical analysis, one finds for 
Msusy £ 2 TeV that m h < m^ ax = m^ ax (tan/5 > 1), where 



m 
m 



h 

max 



122 GeV, if top-squark mixing is minimal, 
135 GeV, if top-squark mixing is maximal. 



(45) 



In practice, parameters leading to maximal mixing are not expected in typical models of supersymmetry 
breaking. Thus, in general, the upper bound on the lightest Higgs boson mass is expected to be 



somewhere between the two extreme limits quoted above. Cross-checks among various programs ||157 
and rough estimates of higher order corrections not yet computed suggest that the results for Higgs 



masses should be accurate to within about 2 to 3 GeV over the parameter ranges displayed in figs. 1^-14. 



In fig. 14, we exhibit the masses of the CP-even neutral and the charged Higgs masses as a function 
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of rriA. The squared-masses of the lighter and heavier neutral CP-even Higgs are related byQ 



m 2 H cos 2 (/3 - a) + m\ sin 2 (/3 - a) = [m^ ax (tan/5)] 2 . (46) 

It is interesting to consider the behavior of the CP-even Higgs masses in the large tan (3 regime. For large 
values of tan/3 and for tua/ tan/3 <C m™ ax (tan/3), the off-diagonal elements of the Higgs squared-mass 
matrix M 2 become small compared to the diagonal elements \M\ 2 \ <C M\]_ + M\ 2 , M\ 2 <C M\ X M\ 2 . 
Hence the two CP-even Higgs squared-masses are approximately given by the diagonal elements of M 2 . 
As above, we employ the notation where m™ ax refers to the asymptotic value of rrih at large tan/3 
and rriA (the actual numerical value of m™ ax depends primarily on the assumed values of the third 
generation squark mass and mixing parameters). If rriA > 7«™ ax , then rrih — m™ ax and m# ~ rriA, 



whereas if rriA < m™ ax , then rrih — vha and rrin ~ m™ ax . This behavior can be seen in fig. [14 



3.2.2 MSSM Higgs mass limits after LEP 



No significant evidence for a Higgs signal has been detected at LEP ||161|| . As a result, one can obtain 
bounds on the possible MSSM Higgs parameters. These limits are often displayed in the m^-tan/3 
plane, although there is additional dependence on various MSSM parameters that effect the radiative 
corrections to the Higgs masses as discussed above. In representative scans of the MSSM parameters, 
the LEP Higgs Working Group |6lJ finds that m h > 91.0 GeV and m A > 91.9 GeV at 95% CL. These 



limits actually correspond to the large tan (5 region in which Zh production is suppressed, as shown in 



fig. 15. In this case, the quoted Higgs limits result of the non-observation of hA and HA 

production. As tan/3 is lowered, the limits on rrih and rriA become more stringent. In this regime, 
the hA production is suppressed while the Zh production rate approaches its SM value. Thus, in this 
case, the SM Higgs limit applies {m h > 114 GeV) as shown in fig. |15](a). The precise region of MSSM 
Higgs parameter space that is excluded depends on the values of the MSSM parameters that control the 
Higgs mass radiative corrections. For example, a conservative exclusion limit is obtained in the maximal 
mixing scenario, since in this case the predicted value of rrih as a function of rriA and tan /3 is maximal 
(with respect to changes in the other MSSM parameters). The excluded regions of the MSSM Higgs 



parameter space based on the maximal mixing benchmark scenario of ref. ||156|| , are shown in fig. |15| , 
and correspond to the exclusion of the range 0.5 < tan/3 < 2.4 at the 95% CL. However, the tan/3 
exclusion region can still be significantly reduced (even to the point of allowing all tan/3 values) by, e.g., 
taking Msusy = 2 TeV and m t = 180 GeV (which still lies within the error bars of the experimentally 
measured value), and allowing for the theoretical uncertainty in the prediction of m™ ax ||162| . 



At tree level, eq. ( fi6l ) is a consequence of the sum rule: J2 n rn ^> n 9^>„zz = \ m z9GGGGi first derived in ref. 131], where 



the tf> n are neutral Higgs bosons of a multi-Higgs-doublet model and ogggg is the quartic coupling of neutral Goldstonc 
bosons. (In the MSSM, gcGGG = fff 2 cos 2 2(3/ cos 2 Qw fl.) A general discussion of related tree-level Higgs mass sum 
rules and bounds can be found in ref. |158| . In fact, cq. (Bq) is more general and applies to the radiatively-corrected 



MSSM Higgs sector [159 160 in the approximation where the renormalizcd a is determined as discussed in Section 3.3.1 
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Figure 15: LEP2 contours of the 95% CL exclusion limits for MSSM Higgs sector parameters as a function of tan/3 and 
(a) mh and (b) tua (in GeV), taken from ref. [161]. The contours shown have been obtained for MSSM Higgs parameters 
chosen according to the maximal mixing benchmark of ref. [1156 



No evidence for the charged Higgs boson has yet been found. The LEP Higgs Working Group quotes 
a limit of m#± > 78.6 GeV at 95% CL ||163 |, which holds for a more general non-supersymmetric two- 
Higgs doublet model and assumes only that the H + decays dominantly into t + v t and/or cs. Although 
the MSSM tree-level bound m H ± > can be relaxed somewhat by radiative corrections, the LEP 
bound quoted above provides no useful additional constraints on the MSSM Higgs sector. 



3.2.3 Effect of explicit CP-violation on the radiatively-corrected MSSM Higgs masses 

In the Standard Model, CP-violation is due to the existence of phases in the Yukawa couplings of the 
quarks to the Higgs field, which results in one non-trivial phase in the CKM mixing matrix. In the 
MSSM, there are additional sources of CP-violation, due to phases in the various supersymmetric mass 
parameters. In particular, the gaugino mass parameters (Mi, i = 1, 2, 3), the Higgsino mass parameter, 
/i, the bilinear Higgs squared-mass parameter, vn\ 2 , and the trilinear couplings of the squark and slepton 
fields (/) to the Higgs fields, Af, may carry non-trivial phases. The existence of these CP phases can 
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significantly affect the MSSM Higgs sector through one-loop radiative corrections [ |164|J165| , |166| , |153|| . 

Note that if one sets fi = Mi = Af = mf 2 , then the MSSM Lagrangian possesses two independent 
global U(l) symmetries — a Peccei-Quinn (PQ) symmetry and an R symmetry. (The quantum numbers 
of the MSSM fields with respect to U(1)pq and U(l)p can be found in refs. ||167f| and [ |166|| .) Conse- 



quently, in the MSSM with nonzero values for the above parameters, there are two independent phase 
redefinitions of the fields that can be used to remove two phases from //, Mi, Af and m\ 2 . However, 
certain combinations of these parameters remain invariant under such phase redefinitions. The simplest 
way to determine these combinations is to treat the aforementioned parameters as spurions with quan- 
tum numbers under the U(1)pq and U(l)p symmetries chosen such that the full MSSM Lagrangian is 
invariant. One can then easily check that the phases of the parameter combinations, aig[fiAf( 



m 



12' 



and axg\piMi(mf 2 )*], are indeed invariant under the U(1)pq and U(l)p phase redefinitions of the MSSM 
fields [|167|Ji66|l . Therefore, if one of these two quantities is different from zero (modulo it), one should 
expect new CP-violating effects induced by the production or exchange of supersymmetric particles. 

We have already noted that the tree-level Higgs sector is CP-conserving. This is a consequence of 
the fact that m\ 2 is the only possible complex parameter that appears in the tree-level Higgs potential. 
Thus the phase of m\ 2 can be rotated away by redefining the phases of the complex Higgs doublets 
appearing in the Lagrangian. The same field redefinition implies that one can choose the vacuum 
expectation values of the two Higgs fields to be real and positive. However, at the one loop-level, 
the Higgs potential acquires a dependence on the parameters \iA t and /iMj through loops of third 
generation squarks and weak gauginos, respectively, which induce non-trivial CP-violating effects. The 
most important of these CP- violating effects is the generation of mixing between the neutral CP-odd and 
CP-even Higgs boson states. Therefore, the physical neutral Higgs bosons are no longer CP-eigenstates 
and the CP-odd Higgs boson mass tua is no longer a physical parameter. The charged Higgs mass is 
still physical and can be used as an input for the computation of the neutral Higgs spectrum of the 



theory ||153|| . The Higgs mass spectrum can therefore be quite different from the CP-conserving case. 
For example, a large splitting between the masses of the next-to-lightest and the heaviest neutral Higgs 
bosons is possible if the charged Higgs boson is not too heavy. 

For large values of the charged Higgs mass, the decoupling limit applies, and the properties of the 
lightest neutral Higgs boson state approach those of the SM Higgs boson. That is, for m H ± ^> mw, the 
lightest neutral Higgs boson is approximately a CP-even state, with CP-violating couplings that are 
suppressed by terms of 0{m^/m l H ±) ||168|| . In particular, the upper-bound on the lightest neutral Higgs 



boson mass, which is reached in the decoupling limit, takes the same value as in the CP-conserving 

case 



||166|| . Nevertheless, there still can be significant mixing between the two heavier neutral mass 
eigenstates. Quantitatively, the leading contribution to the squared-mass terms that mix CP-even and 
CP-odd eigenstates, Mg P (in a convention where m\ 2 is real) is of order 



2 3g 2 mj\fiA 



M sp * c ; , -TF2 (arg[Mtl) ■ (47) 



64vr 2 m^M| 
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Under the reasonable assumption that \fiA t \ < 10MJ, it is clear that the mixing effects between the 
lightest neutral Higgs boson and the heavier Higgs states are small if the masses of the heavy Higgs 
bosons are larger than 2m t . In this limit, the two heavier states are highly degenerate in mass, and 
the CP-violating effects may still lead to non-trivial mixing of the two heavier CP-eigenstates. For 
a detailed study of the Higgs mass spectrum and parametric dependence of the Higgs mass radiative 



corrections, see ref. [153 



3.3 Radiatively-Corrected MSSM Higgs couplings 
3.3.1 Renormalization of cos(/3 — a) 

Radiative corrections also significantly modify the tree-level values of the Higgs boson couplings to 
fermion pairs and to vector boson pairs. As discussed in Section jO, the tree-level Higgs couplings 
depend crucially on the value of cos(/3 — a). In first approximation, when radiative corrections of the 
Higgs squared-mass matrix are computed, the diagonalizing angle a is shifted from its tree-level value. 
Thus, one may compute a "radiatively-corrected" value for cos(/5 — a). This provides one important 
source of the radiative corrections of the Higgs couplings. In fig. [TIL we show the effect of radiative 
corrections on the value of cos(/3 — a) as a function of for different values of the squark mixing 
parameters and tan (3. One can then simply insert the radiatively corrected value of a into eqs. . 

-(|33"D to obtain radiatively-improved couplings of Higgs bosons to vector bosons and to 



([21]) and 
fermions. 

The mixing angle a which diagonalizes the mass matrix in eq. 



0) can be expressed as: 



M 2 



SryC n 
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(Tr.M 2 ) 2 - 4detM 2 
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where s r 



: sin a and c c 
(if Ml < 



cos a. Note that if M\ 2 —> 0, then either sin a — > (if M\ x > 



M 2 
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-M 2 9 ). At tree level, -M 2 2 is small for small and/or large tan/3, but it 



cannot vanish. However, radiative corrections to M. 2 



l 12 = —{m 2 A + m^spCp + 5M.\ 2 can be of the same 
order as its tree level value for small values of and large tan/3. Hence, it is possible for the one- 
loop contribution to approximately cancel the tree-level result (with two-loop corrections to Ai 2 2 small 
compared to the corresponding one- loop result). For moderate or large values of tan/3, the vanishing 
of Ai 2 2 [see eq. (|42]) 1 leads to the approximate numerical relation ||152|| : 



m , 



\ix t tan (3 
lOOMs 



{2a t x t -ll) 



15 
16^ 



In 



m 2 



(49) 



where h t , a s and the weak gauge couplings have been replaced by their approximate numerical values 
at the the electroweak scale. For low values of or large values of the squark mixing parameters, a 
cancellation can easily take place. 
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If M.\ 2 — and tan/3 is large (values of tan/5 > 
couplings is easy to understand. In this limit, M.\i - 



Section 3.2.1 



5 are sufficient), the resulting pattern of Higgs 
t m\ and -M| 2 — m h ax 'j as noted at the end of 
Two cases must be treated separately depending on the value of iua- First, if < m™ ax , 
then sin a ~ —1, cos a ~ and sin/? ~ — cos(/3 — a) ~ 1. In this case, the lighter CP-even Higgs boson h 
is roughly aligned along the $° direction and the heavier CP-even Higgs boson H is roughly aligned along 
the $° direction [see eq. (§)]. In particular, the coupling of H to bb and r + r~ is significantly diminished 
(since down-type fermions couple to $°), while the HVV couplings [eq. (|^)] are approximately equal 
to those of the Standard Model [since cos 2 (/3 — a) ~ 1]. Consequently, the branching ratios of H into 
gg, 77, cc, and W + W~ can be greatly enhanced over Standard Model expectatic 
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160 
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170 



1 and the previous considerations 



Second, if tha 3> m™ ax then sin a ~ and sin (3 ~ cos a ~ sin(/5— a) 
for H apply now to h. 

Although it is difficult to have an exact cancellation of the off-diagonal element M.\ 2 -> m m any 
regions of the MSSM parameter space, a significant suppression of M.\ 2 ma y be present. Generically, 
the leading radiative corrections to M\ 2 depend strongly on the sign of the product \iX t (A t — X t for 
large tan (3 and moderate /j,) and on the value of \At\. For the same value of X t , a change in the sign of \x 
can lead to observable variations in the branching ratio for the Higgs boson decay into bottom quarks. 
If al < 11/2, then the absolute value of M.\ 2 tends to be suppressed [enhanced] for values of fiA t < 
[fj,A t > 0], which implies a similar suppression [enhancement] for the coupling of bottom quarks and 
r-leptons to the SM-like Higgs boson. For larger values of \a t \, the suppression [enhancement] occurs 
for the opposite sign of \iA t . 



3.3.2 The decoupling limit revisited 

Radiative corrections can also significantly affect the onset of the decoupling limit. Recall that at tree 
level [see eq. ([18])], | cos(/3 — a) \ -C 1 for ^> m z , in which case the couplings of h are nearly identical 
to those of the SM Higgs boson. Including the effects of 5A4 2 , we use eq. (H) to obtain 

(M 2 n - M\ 2 ) sin 2(3 - 2Mj 2 cos 2(3 



cos(/3 — a) 



2K 



m 



sin(/9 



a 



m 2 z sin 4/3 + (SM 2 U - 5M 2 22 ) sin 2(3 - 25M 2 l2 cos 2(3 



2(m\ 
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and m\ 

cos(/3 

in the limit of m^ ^> m z , where 
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0(m 2 z ), one finds 
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5M 2 n 
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5M 
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2m| cos 2(3 



m 2 z sin 2(3 



(50) 



(51) 



(52) 
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Eq. ( |5TD exhibits the expected decoupling behavior for ^> mz- However, eq. ( |5UD illustrates another 
way in which cos(/3 — a) = can be achieved — simply choose the MSSM parameters (that govern the 
Higgs mass radiative corrections) such that the numerator of eq. (|50| ) vanishes. That is, 

2m| sin 2/3 = 2 5M\ 2 - tan 2(3 (5M\ X - 8M 2 22 ) . (53) 

Note that eq. (|53|) is independent of the value of rriA- For a typical choice of MSSM parameters, 
eq. ([53|) yields a solution at large tan/3. That is, by approximating tan 2(3 ~ — sin 2(3 ~ — 2/ tan/3, one 
can determine the value of (3 at which the decoupling occurs: 

tm ^M^|±j^. (54) 
oy\4 i2 

The explicit expressions for quoted in eq. ( f4~2"D confirm that the assumption of tan (3 ^> 1 used to 

derive this result is a consistent approximation because 5}A\ 2 ls typically small. We conclude that for 
the value of tan/3 specified in eq. (0), cos(/3 — a) = independently of the value of uia- We shall refer 
to this phenomenon as m^- independent decoupling. From eq. (|4"2"|), it follows that explicit solutions to 
eq. ( |53"D depend on ratios of MSSM parameters and are thus insensitive to the overall supersymmetric 
mass scale, modulo a mild logarithmic dependence on Ms/m t . 



3.3.3 Corrections to tree-level Higgs-fermion Yukawa couplings 



We have seen in Section |3.3.1| that Higgs couplings are modified at one loop due to the renormalization 
of the CP-even Higgs mixing angle a. Additional contributions from the one-loop vertex corrections 
to tree-level Higgs couplings must also be considered ||171| , |172| . |173| . |174j| . These corrections are typically 
small and therefore do not alter significantly the pattern of Higgs couplings. However, at large tan/3, 
the corrections to Higgs-fermion Yukawa couplings can be enhanced, and thus require a careful analysis. 

In the supersymmetric limit, bottom quarks only couple to $° and top quarks only couple to 
However, supersymmetry is broken and a small coupling of the bottom quark [top quark] to $° [$°] 
will be generated from the one-loop Yukawa vertex corrections. These results can be summarized by an 
effective Lagrangian that describes the coupling of the Higgs bosons to the third generation quarks:0 



C 



eff 



5h b )b R <f>' l d Q' L 



(ht + 5h t )i R Q l L <f>i 



Ah t t R Q k L $ k d * 



Ah b b R Q k L <f> k u * 



h.c. 



(55) 



Due to weak isospin breaking, one should allow for different radiatively induced couplings to charged and neutral 
Higgs bosons. For example, one should write AhfJjfjbL^* + A/i^fc at in place of AhbbftQ'l&u* , etc. To the extent 
that weak isospin breaking effects are small in the loop diagrams that generate Ahb and Ahb, it follows that Ahj, f=a Ahf, 
(and similarly for the other radiatively generated coefficients), and we may use eq. ( j5o]) as written. 
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implying a modification of the tree-level relations between h t , h b and m t , m b as follows [ |175| , |176| , |172| . |1 



hv f 5h b Ah b tan(3\ _ h b v . 
m fe = -^=cos/5 I 1 + — + — I = -^=cos/5(l + A 6 ), (56) 

m t = sm ^1 + — + — j— j = ^ sm 0(1 + A t ) . (57) 

The dominant contributions to A;, are tan /5-enhanced, with A b ~ (Ah b /h b ) tan/?; for tan/5 ^> 1, 5h b /h b 
provides a small correction to A b . In the same limit, A t ~ 5h t /h t , with the additional contribution of 
(Aht/ht) cot /3 providing a small correction.^] Explicitly, one finds that for Msusy 3> (where Msusy 
represents a typical supersymmetric mass that appears in the loops) and for tan/5 3> 1 ||175| , |i76| , |144j| , 

|^M § 7(M£, Mf 2 , M|) + I (Ml Ml ^) 



tan /5 , (5£ 



A, - -|^W(M|,M|,M|) - ^/(M^M^ 2 ) , (59) 

where a s = g 2 /4ir, Mg is the gluino mass, M~ bi a are the bottom squark masses, and smaller electroweak 
corrections have been ignored. The loop integral I(a 2 , b 2 , c 2 ) is given by 

^ _ a 2 b 2 \n(a 2 /b 2 ) + b 2 c 2 \n(b 2 /c 2 ) + c 2 a 2 ln(c 2 /a 2 ) 
(a,6,C) " ( a 2_ 6 2 )(6 2_ c 2 )(a 2_ c2) > ( 6 °) 

and is of order l/max(a 2 , 6 2 , c 2 ) when at least one of its arguments is large compared to m 2 z . Note 
that the Higgs coupling proportional to Ah b is a manifestation of the broken supersymmetry in the low 
energy theory; hence, A b does not decouple in the limit of large values of the supersymmetry breaking 
masses. Indeed, if all supersymmetry breaking mass parameters (and /i) are scaled by a common factor, 
the correction A b remains constant. 

Similarly to the case of the bottom quark, the relation between m T and the Higgs-tau-lepton Yukawa 
coupling h T is modified: 

m T = -^(l + A T ). (61) 

The correction A T contains a contribution from a tau slepton-neutralino loop (depending on the two 
tau-slepton masses M^ and Mf 2 and the mass parameter of the B component of the neutralino, Mi) 



21) 



Because the one- loop corrections Sh b , Ahb, Sh t and Ah t depend only on Yukawa and gauge couplings and the 



supersymmetric particle masses, they contain no hidden tan (3 enhancements 
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and a tau sneutrino-chargino loop (depending on the tau sneutrino mass M„ T , the mass parameter of 
the W ± component of the chargino, M 2 , and fi). It is given by ||176| ) 144 |: 



^M 1 ^I(M fl ,M f2 , M 1 ) - ^M 2f iI(M, T ,M 2 ,fi) 



tan/3 



(62) 



where a 2 = 5 ,2 /47r and «i = g' 2 /An are the electroweak gauge couplings. Since corrections to h T are 
proportional to cti and a 2 , they are expected to be smaller than the corrections to hb- 

From eq. ([55]) we can obtain the couplings of the physical Higgs bosons to third generation fermions. 
The resulting interaction Lagrangian is of the form: 



£int 

Using eqs. (|56|) and 



[dhqqhqq + QHqqHqq - ig Aq qAq<y b q} + \bg H - fb tH + h.c. 

q=t,b,T 



(63) 
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(64) 
(65) 
(66) 
(67) 
(68) 
(69) 



and the r couplings are obtained from the above equations by replacing m&, A& and 5hb with m T , 
A T and 5h T , respectively. In addition, one must employ the renormalized value of a in the above 
formulae to incorporate the radiative corrections discussed in section 3.3.1. In writing out the Higgs- 
top quark couplings above, we found it convenient to express the results in terms of A t and Ah t /h t , since 
A t ~ Sh t /h t and the corresponding contribution of Ah t /h t is tan/3 suppressed [eq. fl57|)]. Alternatively, 
eqs. (|6"4])-(|6"9"D can be rewritten in a more symmetrical form by using eqs. (|56D and ( |57D to eliminate A^ 
and A t from the numerators of the corresponding expressions. 
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At large tan /3, terms involving A b oc tan (3 [eq. provide the dominant corrections to the neutral 
Higgs couplings to bb. The corrections proportional to Sh b /h b [see eqs. (|55|) and fl56|) and the discussion 
that follows] are never tan /3-enhanced and are therefore numerically unimportant. The sign of A b is 
governed by the sign of Mgfi, since the bottom-squark gluino loop gives the dominant contribution to 
eq. fl58l). Thus, in a convention where M- g > 0, the radiatively corrected coupling g Abb is suppressed 
(enhanced) with respect to its tree level value for \i > (/1 < 0). In contrast, the radiative corrections 
to 

9hbb an d 9mb have a more complicated dependence on the supersymmetric parameters due to the 
dependence on the CP-even mixing angle a. Since a and A& are governed by different combinations of 
the supersymmetry breaking parameters, it is difficult to exhibit in a simple way the behavior of the 
radiatively corrected couplings of the CP-even Higgs bosons to the bottom quarks as a function of the 
MSSM parameters. 

One can check [using eq. fll9|) 1 that in the decoupling limit, g^q = gh SM qq = Tn q /v. Away from 
the decoupling limit, the Higgs couplings to bottom-type fermions can deviate significantly from their 
tree-level values due to enhanced radiative corrections at large tan/3 [where A b ~ (9(1)] • In particular, 
because A b oc tan/3, the leading one-loop radiative correction to g hb i is of 0{m 2 z tan/3/m^), which 
formally decouples only when m\ 3> m z tan (3. This behavior is called delayed decoupling in ref. [174 



In addition, there are regions of MSSM parameter space in which there is a strong suppression of the 
Higgs coupling to bb (or t + t~) as compared to its tree-level value. As a result, there can be significant 
corrections to the tree-level relation g hb n 1 g hr + T - = m b /m T ||169|| . In particular, in some parameter 
regimes, the r + r~ decay mode can be the dominant h decay channel, a result that would be fatal to 
certain Higgs search strategies which assume that h — > bb is the dominant decay mode. 

Assuming that weak isospin breaking effects in the loop corrections to the charged Higgs fermion 
Yukawa couplings are small (see footnote [19]), then g H - t i [defined in eq. (|63])] is given by 
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(70) 



with a similar form for gn-u T r+ with the replacements noted below eq. 



3.3.4 Effects of explicit CP-violation 

In Section |3.2.3| , we noted the possibility of mixing between the CP-even and CP-odd eigenstates due 
to CP-violating effects that enter via the one-loop radiative corrections. In this case, the neutral scalar 
mass eigenstates, denoted by Hi {i = 1,2,3), are determined by diagonalizing a 3 x 3 squared-mass 
matrix. Thus, one can no longer parameterize the various Higgs couplings in terms of the CP-even 
Higgs mixing angle a. It is convenient to work in a convention where the two vacuum expectation 
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org (A t ) = arg (A b ) [ deg ] arg (A) = arg (A b ) [ deg ] 



Figure 16: (a) Lightest and next-to-lightest neutral Higgs masses and (b) relative couplings (normalized to the SM) of 
the three neutral Higgs bosons to the Z (or W) as a function of the phase of At for the indicated choices of the MSSM 
parameters. Solid [dashed] lines are for arg(Afg) = 0° [90°]. Taken from ref. |153| . 



values are real and positive (by absorbing any potential phases into the definition of the Higgs field) so 
that tan/3 = v u /va as before. Then, eqs. (^) and fl8|) are replaced by 

H t = (v / 2Re$°-^)0 li + (v^Re$° -v u )0 2i + (im $° sin /3 + Im $° cos /?) 3i , (71) 



where O is a 3 x 3 real orthogonal matrix. 

In the CP- violating case, vector boson pairs VV (V = W or Z) couple to all three neutral Higgs 



mass eigenstates, Hi, with ||166 



gHiVV = On cos (3 + 2i sin f3 . (72) 

Fig. [16| shows the dependence of the Higgs masses and the HiZZ squared-couplings on the phase of A t 
for a particular choice of MSSM parameters [as indicated in fig. |T^(a)]. Clearly, these couplings can 
depend sensitively on the phases of the complex supersymmetry-breaking parameters that generate the 
mixing of the CP-even and CP-odd scalar eigenstates through one-loop radiative effects. 

The couplings of V to a pair of neutral Higgs bosons are given by gUiHjzipHi — Ph,), where the 
momenta pn t and pn point into the vertex, gH t H z is antisymmetric under the interchange of Hi and 
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Hj, and ggg ] 



dHiHjZ 



<J 



2 COS 

Using the orthogonality of O (and det O 



[(OwOy - OuOy) sin/3 - (0 3i 2 , - 2i 3i ) cos/3] . 

1), it is easy to derive the relation [116611 :^1 



(73) 



The sum rules of eqs. 



gHiHjZ - Tj^r e ijk g Hk w ■ 
)-(|24T) are then easily extended: 

E2 2 2 

9HiVV ~ 9v m V ' 
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(75) 



which when combined with eq. 
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Ml' 



The couplings of a neutral Higgs boson to H W + are given by gHiH-w+iPHi — Ph-)i where the 



momenta pjj- and point into the vertex, and ||166 

9HiH-W+ — [gHiH+W-]* 



\g [On sin (3 — 2 i cos j3 — iO 



3/ 



(7S 



Another consequence of the CP-violating effects in the scalar sector is that all neutral Higgs scalars 
can couple to both scalar and pseudoscalar fermion bilinear densities (ipip and Vh^, respectively). The 
couplings of the mass eigenstate Hi to fermions depend on the loop-corrected fermion Yukawa couplings, 
hb,ti 5hb,t, Ahb,t, and on tan/3 and the Oji. It is convenient to adjust the phases of the fields so that the 
quantities hb + Shb + Ahbta.n/3 and ht + 5ht + AhtCot/3 are both real and positive [i.e., the physical 
fermion masses are still given by eqs. ( |56"D and (5^)]. The resulting expressions are a straightforward 
generalization of those presented above for the CP-conserving case |153|| : 
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Hff 



i=l 
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b {g%bb + WHibbls) h + -^ t [dHitt + ^Jr.ttTs) * 



(79) 



21 One can easily check that eqs. (20) and (21) are recovered in the CP-conserving limit, where (J7i, #2, H3) = (h, H, A), 
O22 = —On = sin a, O12 = O21 = cos a, O33 = 1, and all other elements Oji vanish. 
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(80) 



(82) 



(83) 



where the Higgs scalar couplings are normalized with respect to the corresponding SM values. 

For large values of the charged Higgs boson mass and for heavy supersymmetric particles, the ex- 
pressions of the lightest neutral Higgs boson coupling to fermions reduce to those of the (CP-conserving) 
SM Higgs boson, as expected for the decoupling limit. In contrast, the two heavy neutral Higgs bosons 
are still admixtures of CP-even and CP-odd eigenstates; hence, CP-violating effects are still present in 
the heavy neutral Higgs sector. However, due to the high degeneracy in mass of the heavy scalar sector 
(especially in the decoupling limit), CP- violating effects may be difficult to observe without precision 
measurements of the heavy neutral Higgs properties. 

The couplings of the charged Higgs bosons to fermions are of the form £ int = b g H - tb tH~ + h.c, 
with 



9H-fb - tan/3 P L 

V v cos p J 



v sm p J 



34) 



One can check that for real Ah b and Ah t , this result is equivalent to eq. ( [TOD [with the same caveats noted 
in footnote |I9||. An explicit computation of the CP- violating H~tb vertex and its phenomenological 



implications can be found in ref. 178 . 
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3.4 MSSM Higgs Boson Decay Modes 



In the MSSM, we must consider the decay properties of three neutral Higgs bosons and one charged 
Higgs pair.0 In the region of parameter space where 3> vrtz and the masses of supersymmetric 
particles are large, the decoupling limit applies, and we find that the properties of h are indistinguishable 
from the SM Higgs boson. If supersymmetric particles are light, then the decoupling limit does not 
strictly apply even in the limit of rriA 3> mz- In particular, the h branching ratios are modified, if the 
decays of h into supersymmetric particles are kinematically allowed. In addition, if light superpartners 
exist that can couple to photons and/or gluons, then the one-loop gg and 77 decay rates would also 
deviate from the corresponding Standard Model Higgs decay rates due to the extra contribution of the 
light superpartners appearing in the loops. In both cases, the heavier Higgs states, H, A and H , are 
roughly mass degenerate, and their decay branching ratios depend crucially on tan (3 as shown below. 

For values of rriA ~ 0(mz), all Higgs boson states lie below 200 GeV in mass. In this parameter 
regime, there is a significant area of the parameter space in which none of the neutral Higgs boson decay 
properties approximates that of the SM Higgs boson. For tan /3 > 1, the resulting Higgs phenomenology 



shows marked differences from that of the SM Higgs boson ||179| . In particular, radiative corrections 
can significantly modify the bb and/or the r + r~ decay rates with respect to those of the SM Higgs 
boson, as noted in Section 3.3.3| . Additionally, the Higgs bosons can decay into new channels, either 



containing lighter Higgs bosons or supersymmetric particles. In the following, the decays of the neutral 
Higgs bosons h, H and A and the decays of charged Higgs bosons are discussed with particular emphasis 
on differences from Standard Model expectations. In the following discussion, we exhibit results for 
tan j3 = 3 and 30 to illustrate the difference between "low" and "high" tan (3. The results shown below 
include the effects of the dominant radiative corrections, which affect both the masses and the couplings 



of the Higgs sector as described in Sections 3.2 and 3.3 



In order to display results for Higgs branching ratios, we must choose a set of MSSM parameters. 
We fix tan f3 (for two representative choices) and vary mj^ from its LEP experimental lower bound of 
90 GeV up to 1 TeV. In addition, the gluino and MSSM squark mass parameters have been chosen to be 
Msusy = Mg = Mq = Mu = Mb = 1 TeV, the squark mixing parameter X t = A t —fi cot (3 = 2.4Msusy, 
and the gaugino mass matrix parameters, /1 = M2 — 2M\ = 1 TeV. This differs somewhat from the 



maximal mixing benchmark scenario of ref. [ 156 1 . Nevertheless, the value of is still close to maximal 
(for fixed and tan/3), so we will continue to loosely refer to the above choice of MSSM parameters 
as a maximal mixing scenario. Our motivation for choosing the gaugino mass parameters large is to 
avoid possible supersymmetric decay modes for the Higgs bosons for Higgs masses below 1 TeV. We 
shall briefly comment on possible supersymmetric decay modes at the end of this section. 

The branching ratios for h and H as a function of their masses are shown in fig. |T7]. As varies 
from 90 GeV to 1 TeV, with the MSSM parameters as specified above, 135 GeV < tuh < 1 TeV when 



22 Unless otherwise noted, we shall neglect CP-violating effects (e.g., by assuming that CP-violating effects induced by 
radiative corrections are small). 
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Figure 17: Branching ratios of the MSSM Higgs bosons h and H, with tan/? = 3 and 30, respectively. Final states labeled 
above include the possibility of one off-shell final state particle below the corresponding two-particle decay threshold. The 
above plots were made under the assumption that the average top and bottom squark masses are f TeV and top-squark 
mixing is maximal. In this case, m™ ax ~ If 5 GeV (f25.9 GeV) for tan/3 = 3 (30), corresponding to the limit of large 
m,A, is indicated by the vertical line in the two left-side plots. The range of run shown corresponds to varying tua 
between 90 GeV and f TeV, while m h > fOO GeV corresponds to m A > 139 GeV (104 GeV) for tan/3 = 3 (30). Other 
supersymmetric parameters have been chosen such that there are no supersymmetric particle decay modes in the Higgs 
mass ranges shown above. 
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Figure 18: (a) Branching ratios of the MSSM Higgs bosons h and H, with tan/3 = 30. Here, we zoom in on the Higgs 
mass regime within ±3 GeV of m™ ax = 125.9 GcV of fig. [l?] in order to get a clearer picture of the various decay modes. 
The range of tuh shown corresponds to 90 GeV < uia < 130 GeV, whereas the range of mu shown corresponds to 
128 GcV < m,A < 1 TcV. (b) Total widths of the MSSM Higgs boson as a function of the corresponding Higgs mass for 
tan j3 = 3 and 30, with the same parameter assumptions employed in fig. [l?]. 



tan/3 = 3 and 126.1 GeV < mjj < 1 TeV when tan/3 = 30. In contrast, most of the variation in 
occurs for values of a few GeV below m™ ax . Thus, we also exhibit in fig. |T8|(a) the branching ratios 
for h and H for tan /3 = 30 and Higgs mass values of m™ ax ± 3 GeV. This reveals a detailed pattern of 
branching ratios that is not easily visible in fig. |17]. The branching ratios for A and H + as a function 
of their masses are shown in fig. [T^. 

The total Higgs decay widths as a function of the corresponding Higgs mass are shown in fig. |l8|(b) 
for the two cases of tan/3 = 3 and 30 (and the other relevant MSSM parameters as described above). 
Note that for large values of the Higgs mass, the corresponding widths are considerably smaller than 
that of the SM Higgs boson. This is due to the suppressed HVV couplings at large Higgs mass and to 
the absence of tree-level AVV and H + W~ Z couplings. One can also check that in the decoupling limit 
{rriA 3> mz), the total width of h coincides with that of /ism- This is illustrated by replotting the h (and 
H) widths on the same plot as the h$M width [see fig. f|(b)]. In particular, note that the dashed and 
dot-dashed h contours in fig. ||(b) approach the /ism contour as reaches its maximal value, (which 
corresponds to the limit of large tua at fixed tan/3). It is interesting to note that in the opposite limit 
of small rriA (especially at large tan/3), cos(/3 — a) — > 1 and it is H that assumes many of the properties 
of hsM- However, there can still be deviations in the Hbb coupling from the corresponding Standard 
Model value at large tan/3, as noted below eq. (^). This explains why the H contours in fig. |](b) 



50 





Figure 19: Branching ratios of the MSSM Higgs bosons A and H + , with tan/3 = 3 and 30, respectively. Final 
states labeled above include the possibility of one off-shell final state particle below the corresponding two-particle decay 
threshold. The above plots were made under the assumption that the average top and bottom squark masses are 1 TeV 
and top-squark mixing is maximal. The range of m H ± shown corresponds to varying mriA between 90 GeV and 1 TeV. 
Other supersymmetric parameters have been chosen such that there are no supersymmetric particle decay modes in the 
Higgs mass ranges shown above. 
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do not quite coincide with the result of the hsM contour as m# approaches its lower limit (with the 
discrepancy between the H and h^u contours more pronounced at large tan/3). 

The branching ratios and widths in figs. |TT|^T§| have been computed using a modified version of 
the HDECAY program ||180|| that incorporates the leading radiative corrections to the Higgs couplings 
discussed in Section The decay modes h,H,A — > bb, t + t~ dominate the neutral Higgs decay 

modes when tan f3 is large for all values of the Higgs masses. For small tan f3, these modes are significant 
for neutral Higgs masses below 2m t (although there are other competing modes in this mass range), 
whereas the ti decay mode dominates above the ti decay threshold. In contrast to the SM Higgs 
boson, the vector boson decay modes of H are strongly suppressed at large mu due to the suppressed 
HVV couplings in the decoupling limit. For the charged Higgs boson, H + — > t + v dominates below tb 
threshold, while H + — ► tb dominates for large values of rriH±- Note that final states labeled in figs. |17| 
and [19] include the possibility of one off-shell final state particle below the corresponding two-particle 
decay threshold [|181|| . For example, for m#± < m t + mj, the tb contour shown in fig. [19| actually 



corresponds to an off-shell t quark that decays to bW + . That is, in this mass region, the tb contour 
corresponds to the branching ratio for the three-body decay H + — > W + bb. This decay mode can be 
especially significant at moderate values of tan f3 due to the large Higgs-top Yukawa coupling. 

As in the Standard Model case, the partial decay widths of the neutral Higgs bosons into bb and cc 
are reduced by about 50-75% when QCD corrections are included (e.g., by employing running quark 
masses in the decay width formulae), whereas the QCD corrections are less significant for Higgs decays 
into ti |29fl . The effects of the QCD radiative corrections on the charged Higgs branching ratios ||182|| are 



significant in the region of tan (3 where the cs and r + u T decay modes are competitive or for large values 
of tan/3 for the decay mode H + — > tb (and for H + — > W + bb below tb threshold [|183|| ). Additional 



super symmetric radiative corrections discussed in Section |3.3| can also significantly affect the Higgs 



boson partial widths. Some of these corrections can be absorbed into the effective mixing angle a [173] 
as shown in Section |3.3.1| . As a consequence of this universal correction, the coupling of h to bb and 
t + t~ can be suppressed for small tua and large tan j3. For the h —>■ bb decay width, the supersymmetric 
corrections | |1 52|]1 71| . |1 72| . |1 73| . |1 74( | proportional to the strong coupling constant a s and the Higgs-top 



quark Yukawa coupling h t can be significant for large values of \i and tan (3. As shown in Section |3.3.3| , 
this effect can be interpreted as a correction to the tree- level relation between m;, and lib- 

In addition to the decay modes of the neutral Higgs bosons into fermion and gauge boson final 
states, there exist new Higgs decay channels that involve scalars of the extended Higgs sector and su- 
persymmetric final states. The unambiguous observation of these modes (as well as any decay mode of 
a charged scalar) would clearly constitute direct evidence of new physics beyond the Standard Model. 
Higgs decays into charginos, neutralinos and third-generation squarks and sleptons can become impor- 



23 For the maximal mixing choice of MSSM parameters used in figs. 17-19[ we find A& ~ 0.55 for tan/3 = 30. For large 
rriA and tan /3, the partial widths of H, A — * bb and H + — > tb (and likewise the corresponding total widths) are suppressed 
by a factor of about (1 + with respect to the corresponding tree-level results. 
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tant, once they are kinematically allowed ||184||. One interesting possibility is a significant branching 



ratio for the decay of a neutral Higgs boson to the invisible mode X1X1 (where Xi is the lightest super- 
symmetric particle). In such a scenario, the discovery of this neutral Higgs boson would be difficult at a 
hadron collider ||185|| . In contrast, at lepton colliders, methods exist for detecting an invisibly decaying 
Higgs boson by observing a peak in the missing mass recoiling against the produced Z. 

3.5 MSSM Higgs Boson Production at Hadron Colliders 
3.5.1 Cross-sections at hadron colliders 

The production mechanisms for the SM Higgs boson at hadron colliders can also be relevant for the 
production of the MSSM neutral Higgs bosons. However, we must take into account the possibility 
of enhanced or suppressed couplings (with respect to those of the Standard Model). For example, 
the HVV couplings are very suppressed in the decoupling limit, and tree-level AVV couplings are 
completely absent. On the other hand, at large tan /3, typically two of the three neutral Higgs couplings 
to bottom-type quarks are enhanced. These effects can significantly modify the neutral Higgs production 
cross-sections. New production mechanisms must be considered for charged Higgs production. 

As in the case of Higgs branching ratios, the predicted cross-sections are sensitive to the MSSM 
Higgs parameters. Again, we consider two representative values of tan/?: a low value of tan/5 = 3 and 
a high value of tan/5 = 30. We then vary ttia, evaluate the other Higgs masses, and compute each 
Higgs cross-section as a function of the corresponding Higgs mass. The Higgs masses and cross-sections 
depend on other MSSM parameters through radiative corrections. As in Section |3.4| , we work in a 
maximal squark mixing scenario in which the value of for a fixed choice of tan (3 and is maximal. 
In addition, because the squark masses are assumed to be heavy (of order 1 TeV), potential supersym- 
metric contributions to the one-loop Higgs-gluon-gluon vertex (due to squark loops) are suppressed. 



Cross-sections for neutral MSSM Higgs production at the Tevatron and the LHC are shown in figs. |20 
and |21] respectively. The dominant Higgs production mechanism over much of the MSSM parameter 
space is gluon-gluon fusion, which is mediated by heavy top and bottom quark triangle loops and the 



corresponding supersymmetric partners [ 186 , 187| , 188|| . The gluon-gluon fusion results shown in figs. 20 



and |2T| include NLO QCD corrections pqfj 

The cross-sections for the production of the neutral CP-even Higgs bosons (0 = h or H) via gauge 
boson fusion V*V* — > <fi (V = W or Z) |45| and via the process qq — > V* — > V<p P2"fl , including first- 



order QCD corrections J 25 ] are also exhibited in figs. and |2l]. Recall that the CP-even scalar <j) has 
SM-like couplings to the vector bosons in two cases: (i) in the decoupling regime for the lightest Higgs 
boson, where </> = h and (ii) for large tan/5 and low m^, where <fi = H. In either case, the SM-like 
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A recent computation of the NNLO QCD corrections to A production via gluon fusion exhibits a 20-30% increase 



over the corresponding NLO cross-section 190 



3 The supersymmetric-QCD corrections due to the exchange of virtual squarks and gluinos are known to be small [188 
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Figure 20: Neutral MSSM Higgs production cross-sections at the Tevatron [y/s = 2 TeV] for gluon fusion gg — > <fi, 
vector-boson fusion qq — > qqV*V* — > qqh, qqH, vector-boson bremsstrahlung qq — > V* — > hV/HV and the associated 
production gg, qq — > (f>bb/ ' <pti including all known QCD corrections, where cf> = h, H or A |3^ , p4| ]. As in fig. |5|, in the vector 
boson fusion process, qq refers to both ud and qq scattering. The four panels exhibited above show (a) h, H production 
for tan/3 = 3, (b) A production for tan/3 = 3, (c) h, H production for tan/3 = 30, (d) A production for tan/3 = 30. 
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Figure 21: Neutral MSSM Higgs production cross-sections at the LHC [y/s = 14 TeV] for gluon fusion gg — > (j>, vector- 
boson fusion qq — > qqV*V* — > qqh, qqH, vector-boson bremsstrahlung qq — > V™ — » hV/HV and the associated production 
99><l<j — > 4>bb/<f)ti including all known QCD corrections, where = fo, i? or A The four panels exhibited above 

show the cross-section in pb vs. the Higgs mass, ranging from 90 GeV to 1 TeV, for (a) h, H production for tan/3 = 3, 
(b) A production for tan/3 = 3, (c) h, H production for tan/3 = 30, (d) A production for tan/3 = 30. For comparison, 
the cross-section for gluon-gluon fusion to a SM Higgs boson is also shown. 
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Higgs scalar, 0, has a mass less than or about equal to 130 GeV, and the corresponding cross-sections 
for V*V* — > and qq — > V* — * V<p are phenomenologically relevant. The other (non-SM-like) CP-even 
scalar has suppressed couplings to VV, and the corresponding cross-sections are generally too small to 
be observed. 

Higgs boson radiation off bottom quarks becomes important for large tan (3 in the MSSM, where the 
Higgs coupling to bottom-type fermions is enhanced. Thus, the theoretical predictions, including full 
NLO computations, are crucial for realistic simulations of the MSSM Higgs signals in these channels.^ 
Moreover, as discussed in Section [3.3.3| , vertex corrections to the bbcf) coupling play a very important 
role in enhancing or suppressing (depending on the MSSM parameters) these production cross-sections 
at large tan (3 W^WMF^ 



We now turn to charged Higgs production. If m#± < m t — mj, then H ± can be produced in the 
decay of the top quark via t — > bH + (and t — > bH~) |192| |. The t — > bH + decay mode can be competitive 
with the dominant Standard Model decay mode, t — > bW + , depending on the value of tan/5, as shown 
in fig. EWa) for m#± = 120 GeV. This figure, taken from ref. | |193j ] illustrates the effects of including 
one-loop radiative corrections. The curved labeled BRqcd, which incorporates the one-loop QCD 
corrections (first computed in ref. [ |194j| ), is applicable to a more general (non-supersymmetric) Type-II 
two-Higgs doublet model [based on the tree- level Higgs- fermion couplings of eqs. (^BJ)-(^)]. Note that 
the supersymmetric corrections can be particularly significant at large tan f3 due to the effect of A& [see 
eq. ( [70^) ] depending on the choice of MSSM parameters. A full one-loop calculation of T(t — > H + b) 
in the MSSM including all sources of large Yukawa couplings can be found in refs. |195|| and ||193|| . A 



treatment including resummation of the leading QCD quantum effects and the dominant contributions 



from loop effects arising from supersymmetric particle exchange can be found in ref. [[196| . 

For m H ± < m t — nib, the total cross-section for charged Higgs production (in the narrow- width 
approximation) is then given by:Q 

a(pp -+H ± + X) = (l- [BR(t -> bW + )f) a(pp -+it + X). (85) 



With a(pp -> tt) ~ 7 pb at y/s = 2 TeV and a{pp -> tt) ~ 1 nb at y/s = 14 TeV PI , roughly 1400 tt 
pairs per detector will be produced per year in Run 2a of the Tevatron (assuming a yearly luminosity 
of 2 fb _1 ), while about 10 7 -10 8 tt pairs will be produced at the LHC (assuming a yearly luminosity of 
10-100 fb _1 ). Folding in the top quark branching ratio, it is a simple matter to compute the inclusive 



2d 



As mentioned in the corresponding discussion for qq, gg — > bbh^M [see Section 2 .2. 1| , one also needs to evaluate 



gb — > b<j) and bb — > <j) with suitable subtraction of the logarithms due to quasi-on-shell quark exchange (to avoid double 
counting) in order to obtain the total inclusive cross-section for <fi producti on. F or example, the tan/3 enhancement can 



lead to copious s-channcl production of Higgs bosons via 6-quark fusion []36 191 ] . 

27 Note that if one evaluates (j(pp — > H + tb + X) in the region of m H ± < m t — m&, one obtains the single charged 
Higgs inclusive cross-section, a(pp — ► H + + X) = BR(i — > bH + )a(pp — » tt + X), rather than full charged Higgs inclusive 
cross-section of eq. (j|q). The latter is not quite a factor of two larger than the former since X can contain a charged 
Higgs boson; one must subtract off [BR(< — ► bH + )] 2 a(pp — -> tt + X) to avoid double-counting. 
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Figure 22: (a) Branching ratio for t — > bH + in the MSSM as a function of tan/3 for m#± = 120 GcV. The three 
curves shown are the results of a computation that (i) is at tree-level; (ii) includes one-loop QCD correction s; an d (iii) 
incorporates both one-loop QCD, electroweak and the effects of MSSM particle exchange (taken from rcf. |193| ); the 
parameters chosen in (iii) correspond to a rather light supersymmetric spectrum: M g = 300 GeV, = 100 GeV, 
M~ bi = 150 GeV, A t = A b = 300 GeV, M 2 = 150 GeV, and M. & = M v = 200 GeV. Curves (i) and (ii) are also applicable 
to a Model-II two-Higgs doublet model without supersymmetry. (b) The charged Higgs production cross-section at the 
LHC near the threshold for t — > bH + for tan/3 = 3 (taken from ref. |199| ). 



charged Higgs cross-section. For values of near m t , the width effects are important. In addition, 
the full 2 — > 3 processes pp — > H + ib+X and pp — > H~tb+X must be considered. In this case eq. (B5|) no 
longer provides an accurate estimate of the charged Higgs cross-section ||198|| , as illustrated in fig. f2^(b) 
(taken from ref. [ 199|| ). The results of fig. |22|(a) imply that for m^± < m t — m^, the discovery of the 



charged Higgs boson at the Tevatron and/or LHC (given sufficient luminosity) is possible if tan/5 ^> 1 
or tan/5 < 1 (the latter is theoretically disfavored). The precise bound on tan/3 (as a function of mjj±) 
depends somewhat on the details of the other MSSM Higgs parameters. 

If rriH± > m t — m^, then charged Higgs boson production occurs mainly through radiation off a 
third generation quark. Single charged Higgs associated production proceeds via the 2 — > 3 partonic 
processes gg, qq — > tbH~ (and the charge conjugate final state). As in the case of bbhsM production, large 
logarithms ]xi(m 2 H ± / m\) arise for m#± ^> m& due to quasi-on-shell t-channel quark exchanges, which 
can be resummed by absorbing them into the 6-quark parton densities. Thus, the proper procedure 
for computing the charged Higgs production cross-section is to add the cross-sections for gb — > tH~ 
and gg — > tbH~ and subtract out the large logarithms accordingly from the calculation of the 2^3 
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Figure 23: The leading-order production cross-sections for charged Higgs production at (a) the Tevatron (pp — > tbH~ + 
X) and (b) the LHC (pp — > tbH~ + X) are shown as a function of ra#± for three values of tan/3 = 2, 10 and 50. The 
cross-sections are obtained by adding the contribution of the 2^2 processes, gb — > tH~ , to those of the 2^3 processes, 
gg — > tbH~ and qq — > tbH~ (suitably subtracted to avoid double counting). Renormalization and factorization scales 
have been both set to m t + m,H±. These results are taken from ref. 201]. 



process ||48| , |200|| . This procedure avoids double-counting of the large logarithms at 0(a s ), and correctly 
resums the leading logs to all orders. In particular, the contribution to the total cross-section coming 
from the kinematical region of the gluon-initiated 2 — * 3 process in which one of the two gluons splits 
into a pair of 6-quarks (one of which is collinear with the initial proton or antiproton), is incorporated 
into the 6-quark parton density. A cruder calculation would omit the contribution of the 2^2 process 
and simply include the results of the unsubtracted 2^3 process. The latter procedure would miss 
the resummed leading logs that are incorporated into the 6-quark density. However, the numerical 
difference between the two procedures is significant only for m#± ^> m t . 



The single inclusive charged Higgs cross-sections at the Tevatron and LHC are exhibited in fig. 23 



as a function of the charged Higgs mass, for tan (3 — 2, 10 and 50. Note that the cross- sect ions shown 
include the region of charged Higgs mass below m#± = m t — corresponding to the case discussed 
above where the charged Higgs cross-section is dominated by ti production followed by t — > bH~ . These 



results are based on the calculations of ref. ||201|| and include the contributions of the 2^2 process 
and suitably subtracted 2^3 process as described above. Similar results have also been obtained in 
ref. [ |202|| . The impact of the leading electroweak and MSSM radiative corrections has been studied in 



ref. [203]. In addition, the NLO QCD corrections to the 2^2 process gb — > H + t have recently been 
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Figure 24: Total cross-section (in fb) for inclusive production of (i) H + H~ +X (207j] (solid line) and (ii) H ± W T +X |205| ] 
(dashed lines) as a function of m H ± for tan/3 = 6 and 30. Curves for (a) pp —> H + H~ + X at the Tevatron and (b) 
pp — > H + H~ + X at the LHC are exhibited. Note that the dependence of process (i) on tan/3 is negligible at the 
Tevatron, while there is some tan /3 dependence at the LHC due to the enhancement of bb — > H + H~ at large tan /3. The 
contribution of bb annihilation to process (ii) dominates over the gg fusion scattering mechanism. AI t = 174.3 GeV and a 
fixed 6-quark pole mass of Mb = 4.7 GeV are used to fix the Higgs-fermion Yukawa coupling. The leading-order CTEQ5L 
parton distribution functions are used. 



evaluated [[204 1. These corrections typically increase the tree-level cross-section by a factor of 1.3 to 1.6, 
depending on the value of the charged Higgs mass and tan /?, with some additional dependence on the 
choice of renormalization and factorization scales. 

Associated production of a charged Higgs boson and a W ± can occur via bb annihilation and gg- 
fusion [205). The contribution of bb annihilation to a{pp — > H^W 1 " + X) [cr(pp — > H ± W T + X)] at the 
Tevatron [LHC] (both charge combinations are included) are shown as function of the charged Higgs 
mass for tan (3 = 6 and 30 in fig. The loop-induced gg fusion contribution is significantly suppressed 
relative to the tree-level bb annihilation if tan/3 > 6, independently of the value of m#± [ 206| . 

Charged Higgs bosons can also be produced in pairs via Drell-Yan qq annihilation. The dominant 
contribution, which arises from uu and dd annihilation into a virtual photon or Z, is independent of 
tan j3. Some tan j3 dependence enters through bb annihilation via t-channel top-quark exchange, although 
this effect is more than one order of magnitude suppressed relative to the dominant contribution at the 
Tevatron. The bb annihilation is more significant at the LHC (at large tan j3 where the H~£b coupling 
is enhanced). The tree-level results for a(pp — > H + H~ + X) at the Tevatron and a(pp — > H + H~ + X) 



at the LHC are shown in fig. 24. These results are obtained |207| with the Higgs-fermion Yukawa 
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coupling based on a fixed 6-quark pole mass of M b 

gg — > H + H~ is typically less important than that of qq annihilation [208 , 209 
in fig 



4.7 GeV. The contribution of the loop-induced 

and is not included 



However, the gg fusion contribution can become significant at large tan (5, with further 
enhancements in some regions of MSSM parameter space in which the squarks (which appear in the 
loop) are light and strongly mixed |p09| . Nevertheless, the inclusive H + H~ cross-section lies below the 
cross-section for single charged Higgs associated production (c.f. figs. ^ and |2~T| ). 

Finally, one can compute the cross-sections for double neutral Higgs production at hadron colliders. 
These include the inclusive production of hh, hH, HH, hA, HA and AA. Cross-sections can be found 
in refs. and | |210| , |211| | (QCD corrections to these cross-sections are evaluated in ref. J55| ). In 

general, the rates for these processes are considerably smaller than for the corresponding single Higgs 
production rates. However, in certain regions of supersymmetric parameter space, squark loops can 
enhance the cross-section for pair production of two CP-even Higgs bosons by as much as two orders 
of magnitude | {210| ]. In some cases, observation of double Higgs production provides some information 
on three-Higgs couplings. For example, for low to moderate values of tan /3, gluon fusion to a virtual 
Higgs boson, which splits into hh, is dominant over bb —> hh. Thus, the overall rate for pp — > hh + X 
would provide a measure of the hhh vertex. 

Additional sources for Higgs boson production can arise from the decay of supersymmetric particles 
into final states containing one or more Higgs boson in the decay chain ||212| . These processes depend 
in detail on the details of the supersymmetric particle spectrum and their couplings. For example, 
the production of h in supersymmetric particle decay followed by the decay h — > bb can yield a signal 
above background at LHC [ 213f| . Processes of this type provide additional channels for possible Higgs 
discovery and precision study, and deserve further analysis. 



3.5.2 Benchmarks for Higgs searches 

In the search for the MSSM Higgs bosons, one must first search for the lightest Higgs scalar, which 
is expected (in almost all cases) to be the neutral CP-even scalar, h. In the decoupling region of the 
MSSM Higgs parameter space (where tua 3> mz), the search techniques already outlined for h$M are 
relevant for h, since the properties of h approximately coincide with those of the SM Higgs boson. The h 
discovery reach can be mapped out as a region of m^-tan j3 parameter space, since these two parameters 
(along with the MSSM parameters that determine the size of the radiative corrections) fix the value of 
m/j. Next, it is critical to identify deviations of the properties of h from those of /ism- Positive evidence 
for such a deviation would signal the existence of additional scalar states of the non-minimal Higgs 
sector. The difficulty of this step depends on how close the model is to the decoupling limit. After the 
discovery of h, the Higgs boson search will focus on the non-minimal Higgs states of the model. 

For values of rriA ~ mz, all the Higgs bosons of the MSSM are of a similar order of magnitude, and 
the properties of h will no longer resemble those of hsM- hi principle, one can then discover multiple 
scalar states in one experiment. Since the two CP-even scalars share the coupling to vector boson pairs 
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[eq. (P^)], one may identify the CP-even scalar whose squared-coupling to VV is larger than 0.5^ gM y V . 
The Tevatron and LHC production cross-sections of this scalar (compared to that of /ism) are reduced 
by no more than 50% (by assumption), while the Higgs branching ratio into bb is similar to that of 
hsM over most of the MSSM parameter space. Thus, the Tevatron and LHC SM Higgs search results 
also apply here modulo minor modifications (which account for the somewhat suppressed production 
cross-section and the effects of supersymmetric corrections to the third generation Yukawa couplings). 

The general MSSM parameter space involves many a priori unknown parameters. In practice, only 
a small subset of these parameters govern the properties of the Higgs sector. Nevertheless, a full scan 
of this reduced subset is still a formidable task. However, a detailed study of a few appropriately 
chosen points of the parameter space can help determine the ultimate MSSM Higgs discovery reach of 
the Tevatron and LHC. It is convenient to choose a set of benchmark MSSM parameters that govern 
the Higgs radiative corrections [ |156|j214| |. These include the supersymmetric Higgs mass parameter fi, 
the third generation squark mixing parameters, A t and Ab, the gluino mass Mg, the diagonal soft- 
supersymmetry-breaking third generation squark squared-masses (which we take for simplicity to be 
degenerate and equal to Msusy), and the top quark mass (which is held fixed at m t = 174.3 GeV). 
The maximal mixing benchmark scenario, is defined as the one in which the squark mixing parameters 
are such that they maximize the value of the lightest CP-even Higgs boson mass for fixed values of 
rriA, tan/3 and M SU sy- Here, we choose X t = A t — /it cot (3 ~ \/Q, Ab = A t , M 2 = — [i = 200 GeV and 
M g = Msusy = 1 TeV [corresponding to m^ ax = 129 GeV]. 

The maximal mixing scenario poses a challenge for Higgs searches, since the predicted Higgs mass 
takes on its maximal value for a given set of MSSM parameters. However, different regions of the MSSM 
Higgs parameter space pose new challenges. For example, regions of parameter space exist in which 
the CP-even neutral Higgs boson with SM-like couplings to the W, Z and t has suppressed couplings 
to bb. The benchmark scenario denoted by "suppressed V<ft — * Vbb production" is an example of this 
behavior. In this case, we take \i = -A t = 1.5 TeV, A b = 0, M 2 = 200 GeV and M g = M SVSY = 1 TeV 
[corresponding to m™ ax = 120 GeV]. The regions of strongly suppressed BR(0 — > bb) correspond to a 
suppressed Hbb coupling at lower m A and a suppressed hbb coupling at larger m^. In particular, the 
suppression for large tan/? extends to relatively large values of ~ 300 GeV, indicating a delay in 
the onset of the decoupling limit. Moreover, in the suppressed V(j) — > Vbb benchmark scenario, all the 
Higgs couplings to bb are generally suppressed, since < <C 1 and sin 2a ~ 0. From the analytic 
formulae, it can be deduced that [iA t < and large values of \A t \, and tan/3 are needed. 

The coverage in the m^-tan/3 plane by different Higgs production and decay channels can vary 
significantly, depending on the choice of MSSM parameters. In the last example in which the CP-even 
Higgs boson with the larger coupling to the W and Z has a strongly suppressed coupling to bottom 
quarks, the Higgs searches at the Tevatron will become more problematical, while the LHC search for 
Higgs production followed by its decay into photons becomes more favorable ||160 . At the same time 



the LHC Higgs discovery reach via vector bosons fusion to Higgs production followed by its decay into 
t + t~ pairs can be significant [215]. 
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3.5.3 MSSM Higgs Boson searches at the Tevatron 



We first consider the Tevatron MSSM Higgs search. Specifically, we make use of the Tevatron hsM search 
techniques, where hsu is replaced by either h or H. If tan/3 ^> 1, a new search mode becomes viable, 
due to the possibility of enhanced couplings of the neutral Higgs boson states to bb |216| , |17Cf| . Thus, we 
also consider the possibility of the bbcf) — > bbbb signature, where <p = h, H, and/or A. If tan/5 is large, 
two of the neutral Higgs boson states, </> = A and h [H] are produced with enhanced rates if < rriz 
[rriA 3> mz], as noted below eq. fl35|). We may combine the results for the various channels to provide 
summary plots of the MSSM Higgs discovery reach of the upgraded Tevatron collider. We consider here 
the results based on a generic MSSM analysis |6"2]] ; see ref. [ 2T7 ] for a similar analysis in the context of a 
variety of models of supersymmetry breaking. In the latter case, the results obtained will be somewhat 
more constraining than the generic analysis, since the supersymmetry-breaking parameters that control 
the radiative corrections to Higgs masses and couplings are no longer arbitrary. 




Figure 25: (a) 95% CL exclusion region and (b) 5a discovery region in the mA~tan/3 plane, for the maximal mixing 
benchmark scenario (see Section 
and gg, qq — ► bb(f> [cf> = h, H, A], <f> — > bb (region in the upper left-hand corner bounded by the solid lines; the two sets 
of lines correspond to CDF and D0 simulations). The region below the solid black line near the bottom of the plot is 
excluded by the absence of e + e~ — > Z<f) events at LEP. Taken from ref. p2j . 



3.5.2) and two different search channels: qq — » V(f> [<j> = h, H], <j> — » bb (shaded regions) 
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Figure 26: The same as fig. |25| but for the suppressed V(f> — ► Vbb production benchmark scenario of Section |3.5.2| . Taken 
from ref. |62|. 



In figs. ^5] and |26], we show the regions of 95% CL Higgs exclusion and 5cr Higgs discovery on the 
m^-tan (3 plane, for two representative MSSM parameter choices, via the search of neutral Higgs bosons 
in the channels: qq — > V(p [<p = h, H], — > bb (shaded regions) and gg, qq — > bbcp [<p = h, H , A], — ► bb 
(region in the upper left-hand corner bounded by the solid lines), for different integrated luminosities 
as indicated by the color coding. The shaded regions presented in these figures reflect the results of 
the SHW simulation of qq — > V<p improved by neutral network techniques |H| . The two sets of lines 
(for a given color) bounding the regions accessible by the bbcf) search correspond to the CDF and D0 
simulations, respectively. The solid black line near the bottom of each plot indicates the lower limit of 
tan j3 (as a function of m^) based on the absence of observed e + e~ 
the importance of the complementarity between the qq — > Vcf) and qq 



1611. Note 



Zcf) events at LEP 
bb(p channels for improving the 
coverage of the MSSM parameter space in the low rriA region in fig. ^6|(a). The results of figs. |25|(a) 
and |26"l(a) demonstrate that 5 fb -1 of integrated luminosity per experiment will be sufficient to cover 
nearly all of the MSSM Higgs parameter space at 95% CL in the benchmark scenarios specified above. 

To assure discovery of a CP-even Higgs boson at the 5a level, the luminosity requirement becomes 
very important. Figs. |25|(b) and ^6|(b) show that a total integrated luminosity of about 20 fb -1 per 
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experiment is necessary in order to assure a significant, although not exhaustive, coverage of the MSSM 
parameter space. In general, we observe that the complementarity between the two channels, qq — * V<p 
and qq — > bb(f>, is less effective in assuring discovery of a Higgs boson as compared with a 95% CL Higgs 
exclusion. This is due to the much higher requirement of total integrated luminosity combined with 
the existence of MSSM parameter regimes which can independently suppress both Higgs production 
channels. Fig. 2t exhibits one of the most difficult regions of MSSM parameter space for Higgs searches 
at the Tevatron collider. Nevertheless, even in this case, a very high luminosity experiment can cover 
a significant fraction of the available MSSM parameter space. 

If explicit CP violation occurs through nonzero phases of the supersymmetry breaking parame- 
ters, then the three neutral Higgs bosons are a combination of CP-even and CP-odd states and the 
phenomenology can become much more complicated. In particular the couplings of the neutral Higgs 
bosons to the W and Z bosons are now shared by the three Higgs bosons and it may well be that the 
lightest Higgs has such a weak coupling to the vector bosons that it would have been missed at LEP and 
will be elusive at the Tevatron. Fig. ^ shows an interesting example where the effects of CP violation 
are such that for CP- violating phases of the parameter A t of about 90°, the lightest Higgs boson cannot 
be detected at the Tevatron even though its mass is below 100 GeV, but the second lightest Higgs has 
SM-like couplings to the W and Z and thus can be detected if sufficient luminosity is provided. 



3.5.4 MSSM Higgs Searches at the LHC 

If no Higgs boson is discovered at the Tevatron, the LHC will cover the remaining unexplored regions of 



the m^-tan (3 plane, as shown in figs. |27] and |28|. That is, in the maximal mixing scenario (and probably 
in most regions of MSSM Higgs parameter space), at least one of the Higgs bosons is guaranteed to 
be discovered at either the Tevatron and/or the LHC. A large fraction of the parameter space can be 
covered by the search for a neutral CP-even Higgs boson by employing the SM Higgs search techniques, 
where the SM Higgs boson is replaced by h or H with the appropriate re-scaling of the couplings. 
Moreover, in some regions of the parameter space, both h and H can be simultaneously observed, and 
additional Higgs search techniques can be employed to discover A, and/or at the LHC. 

A CP-even Higgs boson, <p, can be observed in a number of different decay modes. If > m™ ax , 
then <fi = h is SM-like (near the decoupling limit), whereas at large tan/3 and < m™ ax , <ft = H 
is the SM-like Higgs bosonQ It is possible to observe <fi — > 77 when <fi is produced singly via gg and 
fusion, or when produced in association with W and/or tt. A second decay mode, h — > bb, 
can be observed in <ptt production. Finally, it may be possible to observe </> — > t + t~ when is 



produced via V*V* fusion, where the forward jets are used to help reduce backgrounds [215]. Hence 



by using the complementarity of the various Higgs signatures described above, one can discover over 



28 For vtia < m™ ax and moderate tan j3 values, neither CP-even Higgs boson is SM-like, although both Higgs masses lie 
below about 150 GeV and will appear (albeit with reduced couplings to VV) in the Higgs searches described above. 
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Figure 27: (a) 5<r discovery contours for MSSM Higgs boson detection in various channels in the m^-tan/3 plane, in 



the maximal mixing scenario, assuming an integrated luminosity of L = 30 fb for the CMS detector [ 218 1 . (b) As in 



(a), but for an integrated luminosity of L = 300 fb 1 for the ATLAS detector [ 222 



nearly the entire MSSM parameter space, given sufficient integrated luminosity.^] In order to illustrate 
the complementarity of the 77 and bb decay modes, we exhibit in fig. ^ the regions of MSSM Higgs 
parameter space that can be covered for the two benchmark scenarios of MSSM parameters described in 
Section |3.5.2| . The behavior illustrated in this figure can be understood by noting that the <pbb coupling 
can be significantly suppressed (or enhanced), depending on the impact of the radiative corrections 
discussed in Section p.3| . As a result, the branching ratio for <fi — > 77 is correspondingly larger (or 
smaller), with obvious implications for the <p — > bb and — > 77 searches. 

We next focus on the potential for observing the heavier Higgs states {H^, A and H). A number of 
recent studies |74| . |75| . |2T8| . [2 19| . p20|J22 1|| show that the following modes will be effective in searching for the 
heavier MSSM Higgs bosons. For the heavy neutral Higgs bosons, the most relevant decay signatures 



29 One must still demonstrate that it is possible at the LHC to discover the lightest CP-even Higgs boson, even if its 
branching ratios into bb and/or 77 are significantly suppressed (either due to the effects of radiative corrections or due 
to the existence of a significant branching fraction into invisible modes). Such suppressions can occur in regions of the 
MSSM parameter space not yet considered by the LHC Higgs search simulations. 
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Figure 28: Regions in the m^-tan (3 plane in the maximal mixing scenario in which up to four Higgs boson states of the 
MSSM can be discovered at the LHC with 300 fb _1 of data, based on a simulation that combines data from the ATLAS 
and CMS detectors. Taken from ref. IB7II. 



are: A, h — > t + t~ (where the r is detected either via its leptonic or hadronic decay) and A, h — > 
which yield promising signals if tan j3 is large. The r + r~ channel provides the largest discovery reach 
in the heavy Higgs mass. Other possible neutral Higgs decays: A, H —> ti; H — > ZZ* — > A£; H — > hh 
and A — > Zh are significant in regions of the parameter space that are (nearly) ruled out by the LEP 
Higgs search. For the charged Higgs boson, we must again consider whether can be produced 
in (on-shell) top-quark decays. If this decay is forbidden, the positively charged Higgs boson will be 
produced primarily by gb — > H + t (see Section p.5.1[) . In either case, the observation of the charged 
Higgs boson is possible if tan/3 3> 1 or tan/3 < (9(1) ||221 



For large tan/3, the decays H~ 



► r 1 v 

and tb (if kinematically allowed) provide the most favorable signatures. In particular, the rv decay 
mode, followed by the hadronic decay of the r provides the largest discovery reach for large m#±. The 
ultimate charged Higgs mass reach can depend significantly on the choice of MSSM parameters that 
control the radiative corrections to the Higgs-bottom quark Yukawa coupling [|196|1 [see, e.g., eq. (fT0|)1. 

Putting all of the above results together, it may be possible at the LHC to either exclude the 
entire m^-tan/3 plane (thereby eliminating the MSSM Higgs sector as a viable model), or achieve a 
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Figure 29: Complementarity between the LHC searches for the decay modes — > bb and — ► 77, where <f> = h or H 
corresponds to the CP-even Higgs boson with the larger coupling to VV. Two different choices of MSSM parameters 
are exhibited: (a) the maximal mixing scenario and (b) the suppressed Vb — > Vbb production scenario. In both cases, 
the region corresponding to 5cr discovery of 4> ~ > 77 with 30 fb _1 of data is shaded with parallel horizontal lines. With 
100 fb _1 of data, these regions expand to include areas shaded with diagonal parallel lines with positive slope. The region 
corresponding to 5a discovery of ti<p — > tibb with 30 fb~ x of data is shaded with diagonal parallel lines with negative 
slope. With 100 fb _1 of data, these regions expand to include the entire m^-tan/3 plane excluding the blue cross-hatched 
region. The "unexcluded" region (where no discovery of ticj> — + tibb is possible) occupies a small region at low tan (3 and 
m A in both (a) and (b). In addition, in case (b), the excluded region also includes the two narrow wedge regions at large 
tan (3 and low tua- Taken from ref. |223 



5cr discovery of at least one of the MSSM Higgs bosons, independently of the value of tan/3 and m A . 
For example, Fig. [27] shows what can be achieved by the CMS detector with 30 fb _1 ||218|| and by the 
ATLAS detector with 300 fb^ 1 [ 222 1 , assuming the maximal mixing scenario. Note that over a significant 



fraction of the MSSM Higgs parameter space, at least two Higgs bosons can be observed as shown in 
fig J 



_|. Nevertheless, there is still a sizable wedge-shaped region at moderate values of tan (3 opening 
up from about = 200 GeV to higher values in which the heavier Higgs bosons cannot be discovered 
at the LHC. In this parameter regime, only the lightest CP-even Higgs boson can be discovered, and 
its properties are nearly indistinguishable from those of the SM Higgs boson. Precision measurements 
of Higgs branching ratios and other properties will be required in order to detect deviations from SM 
Higgs predictions and demonstrate the existence of a non-minimal Higgs sector. 
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Finally, we noted at the end of Section |3.5.3| that CP-violating effects in the Higgs sector can 
modify the usual CP-conserving Higgs phenomenology. As a result, the LHC discovery reach of various 
Higgs channels discussed above may be altered in a significant way. It is therefore essential to make 
complementary measurements in as many Higgs channels as possible in order to cover the most general 
MSSM parameter space [|223|| . 



3.6 MSSM Higgs Boson Searches at the LC 



The main production mechanisms for the MSSM Higgs bosons are ||22 1 



(i) e + e~ — > Zh , ZH via Higgs-strahlung , 

(ii) e + e~ — > vvh , vvH via W + W~ fusion , 

(Hi) e + e~ — > hA , HA via s-channel Z exchange , 

(iv) e + e~ — > H + H~ via s-channel 7 , Z exchange . (86) 

As in the SM Higgs search, process (i) followed by Z — > £ + £~ allows the Higgs boson recoiling against 
the Z to be reconstructed, independently of the Higgs decay channel. Thus, h [or H if sin(/5 - a)<l] 
can be discovered at the LC even if it has a large branching fraction into invisible modes.Q 

Processes (i) and (Hi) are complementary to each other as a consequence of unitarity sum rules for 
tree-level Higgs couplings ||131| , |132|| . In particular, eq. ( [24]) implies that both g\zz an d 9\az (0 = h 



or H) cannot simultaneously vanish. If < m™ ax , then all the MSSM Higgs boson states have 
mass below 150 GeV, and can be cleanly reconstructed at the LC (with yfs > 350 GeV) via the four 
production mechanisms listed above [^| . On the other hand, when > 200 GeV, one finds that 
rriA ~ ran ~ rn H ± and gnzz ~ 9haz ~ 0, and the couplings of h are nearly identical to those of 
hsm as a consequence of the decoupling limit. Since nih < 135 GeV, the LC with a center-of-mass 
energy of 300 GeV is more than sufficient to observe the h [via processes (%) and (ii)] and thus cover 
the entire MSSM parameter space with certainty. Moreover, the cross-sections for HZ, Hvv and hA 
are strongly suppressed [since | cos(/3 — a)\ 1]. The cross-sections for HA and H + H~ production 
are unsuppressed if kinematically allowed.^ That is, the heavy Higgs bosons, H, A and H ± can only 
be observed in pair production processes where both Higgs states are heavy (and the minimum y'i 
required is somewhat above 2m a) ■ These features are evident in fig. |3(| which depicts cross-sections 
for Higgs-strahlung [process (i)\ and associated Higgs pair production [processes (Hi) and (iv)] as a 
function of the corresponding Higgs mass for two different choices of y/s and tan (3. The cross-section 

30 The processes e + e~ — > e + e~h, e + e~H (via ZZ fusion) also allow for Higgs detection, independently of its decay 
channel, by reconstruction of the Higgs boson recoiling against the final state e + e~ pair. However, the ZZ fusion rates 
are an order of magnitude smaller than the corresponding W + W~ fusion rates [process (ii)]. 

31 Duc to the p-wave suppression at threshold, the HA and H + H~ cross-sections fall off rapidly as the corresponding 
Higgs masses approach yfs/2. 
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Figure 30: MSSM Higgs boson production rates at the LC for two choices of tan/? = 3 (solid) and 30 (dotted) for 
(a) y/s = 350 GeV as a function of the mass of the produced CP-even neutral Higgs boson (either h or H); and for 



(b) y/s = 800 GeV as a function of m H ± and tua, respectively. Taken from ref. [225 



for Higgs production via W + W fusion [process (ii)] is not shown. The uucj) production cross-section 
is suppressed relative to the corresponding SM cross-section (shown in fig. |9|) by a factor of sin 2 (/3 — a) 
[cos 2 ((3-a)] for = h [<p = H). 

In addition to H + H~ production, there are a number of mechanisms in which the charged Higgs 
boson is singly produced. Charged Higgs bosons can be produced in top decays via t —> b + H + if 
mn± < m t — m&, as discussed previously in Section |3.5.1 . The process e + e~ —>■ W ± H T , which arises 



at one-loop [ |226| , |227| , allows for the possibility of producing a charged Higgs boson with m#± > y/s/2, 
when H + H~ production is kinematically forbidden. With favorable MSSM parameters and moderate 
values of tan/3, more than ten W^H^ events can be produced at the LC for ra#± < 350 GeV with 
yfe = 500 GeV and 500 fb" 1 of data, or for m H ± < 600 GeV with yfs = I TeV and 1 ab" 1 p7| . Other 



single charged Higgs production mechanisms include tbH~ / ibH + production | 101 |, r + vH~ / t~vH + 
production [ |228| ] , and a variety of processes in which if^ is produced in association with a one or two 
other gauge and/or Higgs bosons ||229|| . 



The heavier Higgs states could lie beyond the discovery reach of the LC (\/s < 1 TeV) and the LHC 
[cf. fig. |28|]. In this case, the precision measurements of the h decay branching ratios and couplings 
achievable at the LC are critical for distinguishing between h$M and h of a non-minimal Higgs sector 
with properties close to that of the SM Higgs boson. To illustrate the challenge of probing the decoupling 
limit, suppose that > a/s/2 so that only the light Higgs boson, h, can be observed directly at the 
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Figure 31: Contours of \ 2 f° r Higgs boson decay obscrvablcs for (a) the maximal mixing scenario; and (b) a choice 
of MSSM parameters for which the loop-corrected hbb coupling is suppressed (relative to the corresponding tree-level 
coupling) at large tan/3 and low ttia- These results are based on Higgs partial width measurements anticipated at the LC 
(shown in Table 1) with i/s = 500 GeV and an integrated luminosity of 500 fb _1 . The contours correspond to 68, 90, 95, 



98 and 99% confidence levels (right to left) for the observables g\ hh , 9h TT , an d ff/L„- See rcf. [ 1 1 2 1 for additional details 



LC. In this case, the fractional deviation of the couplings of h relative to those of the SM Higgs boson 
scales as m 2 z /m 2 A . Thus, if precision measurements reveal a non-zero deviation, one could in principle 
derive a constraint {e.g., upper and lower bounds) on the heavy Higgs masses of the model. In the 
MSSM, the constraint is sensitive to the MSSM parameters that control the radiative corrections to 
the Higgs couplings. This is illustrated in fig. BTL where the constraints on m A are derived for two 
different sets of MSSM parameter choices ||112|| . Here, a simulation of a global fit of measured hbb, hrr 
and hgg couplings is made (based on the results of Table 1) and x 2 contours are plotted indicating the 
constraints in the m^-tan f3 plane, assuming a deviation from SM Higgs boson couplings is seen.0 In the 
maximal mixing scenario shown in fig. [H](a), the constraints on are significant and rather insensitive 
to the value of tan/3. Similar results, in which deviations of BK(h — > bb)/BK(h — > t + t~) from the 



32 These results are similar to those of ref. [107], in which rrih = 120 GeV and an integrated LC luminosity of 1000 fb _1 
were assumed. The upper bound on uia as a function of tan/3 was determined in which 68%, 90% and 95% of the 
MSSM parameter space, respectively, yield h branching ratio predictions that differ from those of the SM at th e 95 % CL. 
However, the MSSM parameter regimes in which the sensitivity to ttla are weakest were not identified in ref. [ 107 1 . 
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SM prediction yield limits on the allowed values of rriA as a function of tan/3, have been obtained in 
ref. ||230|| . However in some cases, as shown in fig. |31|(b), a region of tan (3 may yield almost no constraint 
on rriA- This is due to the phenomenon of m^-independent decoupling noted below eq. (|54|), in which 
cos(/3 — a) [which controls the departure from the decoupling limit] vanishes at a particular value of 
tan/3 independently of the value of tua- Thus, one cannot extract a fully model-independent upper 
bound on the value of (beyond what can be deduced if no direct A production is observed at the 
LC). Of course, after supersymmetric particles are discovered, information about the MSSM spectrum 
can be used to obtain a more stringent bound on rriA, using the techniques described above. 

The e + e~ linear collider running in the 77 collider mode presents additional opportunities for the 
study of the MSSM Higgs sector. Resonance production 77 — > H and A can be used to extend the 



reach in Higgs masses beyond the limit set by HA pair production in the e + e mode 1 231, 123, 124 



Typically, one can probe the heavy Higgs masses out to ~ 0.8-y/i (where yfs is the center of mass 
energy of the LC). This expands the MSSM Higgs discovery reach further into regions of the m^-tan/3 
parameter space for which the LHC is not sensitive in general (the so-called "blind wedge" of large 
and moderate values of tan (3 seen in fig. ^8|) . 

As noted above, at least one Higgs boson must be observable at the LC in the MSSM. In non-minimal 
supersymmetric models, additional Higgs bosons appear in the spectrum, and the "no-lose" theorem 
of the MSSM must be reconsidered. For example, in the non-minimal supersymmetric extension of 
the Standard Model (the so-called NMSSM where a Higgs singlet is added to the model ||232|| ), the 
lightest Higgs boson decouples from the Z boson if its wave function is dominated by the Higgs singlet 
component. However, in this case the second lightest neutral CP-even Higgs boson usually plays the 
role of h of the MSSM. That is, the mass of the second lightest Higgs boson is typically below 150 GeV 
with significant couplings to ZZ , so that it can be produced by the Higgs-strahlung process with an 
observable cross-section ||233| , |158|| . If the second lightest Higgs boson also decouples from the Z, then the 
third lightest Higgs boson will play the role of h of the MSSM for which the observation is ensured, and 
so on. Even in bizarre scenarios where all the neutral Higgs boson share equally in the coupling to ZZ 
(with the sum of all squared couplings constrained to equal the square of the h^ZZ coupling [|131|J132|| ), 
the "no-lose" theorem still applies — Higgs production at the LC must be observable |[234|| . In contrast, 
despite significant progress, there is no complete guarantee that at least one Higgs boson of the NMSSM 



must be discovered at the LHC for all choices of the model parameters ||235 



One of the key parameters of the MSSM Higgs sector is the value of the ratio of Higgs vacuum 
expectation values, tan/3. In addition to providing information about the structure of the non-minimal 
Higgs sector, the measurement of this parameter also provides an important check of supersymmetric 
structure, since this parameter also enters the chargino, neutralino and third generation squark mass 
matrices and couplings. Thus, tan/3 can be measured independently using supersymmetric processes 
and compared to the value obtained from studying the Higgs sector. Near the decoupling limit, the 
properties of h are almost indistinguishable from those of /ism ; and thus no information can be extracted 
on the value of tan/3. However, the properties of the heavier Higgs bosons are tan (3- dependent. Far 



71 



from the decoupling limit, all Higgs bosons of the MSSM will be observable at the LC and exhibit 
strong tan /3-dependence in their couplings. Thus, to extract a value of tan j3 from Higgs processes, one 
must observe the effects of the heavier Higgs bosons of the MSSM at the LC. 

The ultimate accuracy of the tan/5 measurement at the LC depends on the value of tan/5. In 



Ref. ||236|| , it is argued that one must use a number of processes, including bbbb final states arising from 
bbH, bbA, and HA production, and tibb final states arising from tbH + , biH~ and H + H~ production. 
One subtlety that arises here is that in certain processes, the determination of tan j3 may be sensitive 
to loop corrections that depend on the values of other supersymmetric parameters. One must settle on 
a consistent definition of tan/3 when loop corrections are included |[237|| . A comprehensive analysis of 
the extraction of tan (3 from collider data, which incorporates loop effects, has not yet been given. 

The study of the properties of the heavier MSSM Higgs bosons (mass, width, branching ratios, 
quantum numbers, etc.) provides a number of additional challenges. For example, in the absence of 
CP-violation, the heavy CP-even and CP-odd Higgs bosons, H and A, are expected to be nearly mass- 
degenerate. Their CP quantum numbers and their separation can be investigated at the same time in 
the 77 collider mode of the LC. If the polarization states of the two incoming linearly-polarized photons 
are parallel [perpendicular] then only the CP-even Higgs boson H [CP-odd Higgs boson A] will be 
produced ||238|| . Thus, the determination of the Higgs boson CP quantum numbers and the separation 
of the two different states can be achieved. In the case of a CP-violating Higgs sector, the observation 
and measurement of Higgs boson properties become much more challenging. The 77 collider can provide 
new opportunities to test the nature of the couplings of the Higgs neutral eigenstates (with indefinite 
CP quantum numbers) to gauge bosons and fermions ||124| , |239|| . 

Finally, once the heavy Higgs spectrum is revealed, one would like to reconstruct the two-Higgs- 
doublet scalar potential ||240|| . This is not likely to be accomplished at a first generation LC, although one 
can make a start if the heavy Higgs masses are not too large. To probe aspects of the Higgs potential one 
must observe multiple Higgs production in order to extract the Higgs self-couplings [ 113| , 114 , 240 , p4l 



Ultimately, such a program would require an LC with very high energy and luminosity such as CLIC. 



4 Conclusions 

The physical origin of electroweak symmetry breaking is not yet known. In all theoretical approaches 
and models, the dynamics of electroweak symmetry breaking must be revealed at the TeV-scale or 
below. This energy scale will be thoroughly explored by hadron colliders, starting with the Tevatron 
and followed later in this decade by the LHC. Even though the various theoretical alternatives can 
only be confirmed or ruled out by future collider experiments, a straightforward interpretation of the 
electroweak precision data suggests that electroweak symmetry breaking dynamics is weakly-coupled, 
and a Higgs boson with mass between 100 and 200 GeV must exist. With the supersymmetric extension 
of the Standard Model, this interpretation opens the route to grand unification of all the fundamental 
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forces, with the eventual incorporation of gravity in particle physics. 

In this review we have summarized the theoretical properties of the Standard Model Higgs boson 
and the Higgs bosons of the MSSM, and surveyed the search strategies for discovering the Higgs boson 
at hadron and lepton colliders. We have assessed the Higgs boson discovery reach of present and future 
colliders, and described methods for measuring the various Higgs boson properties (mass, width, CP 
quantum numbers, branching ratios and coupling strengths). 

The observation of a Higgs boson in the theoretically preferred mass range below 200 GeV may be 
possible at the Tevatron, whereas experiments at the LHC can discover the SM Higgs boson over the 
full Higgs mass range up to 1 TeV. The Tevatron can also extend the LEP search for Higgs bosons of 
the MSSM by either discovering the lightest CP-even MSSM Higgs boson, h (or in some special cases 
discovering additional Higgs scalars of the model), or by further constraining the MSSM Higgs parameter 
space. The LHC is sensitive to nearly the entire MSSM Higgs parameter space, in which either h alone 
can be discovered or multiple Higgs states can be observed. A program of Higgs measurements will be 
initiated at the LHC to measure Higgs partial widths with an accuracy in the range of 10-30%. 

The discovery of the Higgs boson at the Tevatron and/or the LHC is a crucial first step. The 
measurement of Higgs properties at the LHC will begin to test the dynamics of electroweak symmetry 
breaking. However, a high-luminosity e + e~ linear collider, now under development, is needed for a 
systematic program of precision Higgs measurements. For example, depending on the value of the 
Higgs mass, branching ratios and Higgs couplings can be determined in some cases at the level of a few 
percent. In this way, one can extract the properties of the Higgs sector in a comprehensive way, and 
establish (or refute) the existence of scalar sector dynamics as the mechanism responsible for generating 
the masses of the fundamental particles. 
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